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FINAL  REPORT 


GRANT#:  N00014-05-1-0009 

PRINCIPAL  INVESTIGATORS:  Diola  Bagayoko,  Ph.D.,  Southern  University  System 
Distinguished  Professor  of  Physics,  Project  Director,  and  G.  L.  Zhao,  Ph.D.,  Associate  Professor 
of  Physics,  SUBR,  Co-Principal  Investigator 

Program  Officer,  Office  of  Naval  Research  (ONR):  Colin  Wood,  DSc,  Ph.D.,  Code  312, 
Office  of  Naval  Research  (ONR),  Room  9-06,  1  Liberty  Center,  875  N.  Randolph  St.,  Arlington 
VA  22203.  Telephone:  (703)  696-4218. 

INSTITUTION:  Southern  University  and  A  &  M  College  in  Baton  Rouge  (SUBR) 

GRANT  TITLE:  Predictive  Computations  of  Properties  of  Wide-Gap  and  Nano- 
Semiconductors 

AWARD  PERIOD:  October  1,  2004  to  December  31,  2006 

OBJECTIVE:  As  per  the  initial  proposal,  the  aim  of  this  project  is  (a)  to  produce  new  scientific 
knowledge  on  wide-gap  and  nano-semiconductors,  with  specific  materials  including  GaN,  InN, 
AIN,  and  carbon  nanotubes,  (b)  to  develop  analytical  and  computational  techniques  and 
computer  codes  for  wider  applications  in  descriptive  and  predictive  investigations  of  the  above 
materials  and  others,  and  (c)  to  provide  research  training  for  undergraduate  and  graduate 
students,  including  minority  students  at  SUBR. 

APPROACH:  Our  technical  approach  or  method,  as  per  the  proposal,  consists  of  the  following. 
Like  most  previous  calculations  of  electronic  and  related  properties  of  semiconductors,  we 
employ  a  density  functional  theory  (DFT)  potential.  This  potential  is  often  selected,  by  design, 
from  the  local  density  approximation  (LDA)  to  the  general  DFT.  We  also  utilize,  like  previous 
works,  the  linear  combination  of  atomic  orbital  (LCAO)  formalism  in  a  Rayleigh-Ritz  variational 
procedure  to  perform  self-consistent  calculations.  Our  calculations,  to  date,  are  for  zero 
temperature  and  are  non-relativistic. 

The  key  distinction  between  our  method  and  those  of  all  previous  calculations  by  other 
groups  resides  in  our  adherence  to  the  Bagayoko,  Zhao,  and  Williams  (BZW)  method  in  the 
implementation  of  the  LCAO  formalism.  References  clearly  identified  in  the  attached  articles 
fully  describe  this  method.  Essentially,  this  method  resolved  the  long-standing  underestimation, 
by  theory  (by  30  to  50%  or  more),  of  measured  energy  or  band  gaps  of  finite  (atoms,  clusters, 
molecules)  and  infinite  (semiconductors  insulators)  systems.  As  explained  in  the  attached  articles 
and  references  therein,  BZW  first  identified  a  basis  set  and  variational  effect  that  is  responsible 
for  an  unphysical  lowering  of  the  low-lying,  unoccupied  energy  levels  or  bands  of  materials. 
This  unphysical  lowering  is  shown  to  be  a  mathematical  artifact  stemming  from  the  Rayleigh 
theorem.  BZW  subsequently  introduced  a  rigorous  method  of  performing  variational 
calculations  in  such  a  way  that  the  above  unphysical  lowering  of  unoccupied  levels  or  bands  is 
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totally  avoided.  Let  us  underscore  that  the  underestimation  of  energy  and  band  gaps  is  a  trivial 
consequence  of  the  above  unphysical  lowering  of  unoccupied  levels  of  bands. 

The  BZW  method  generally  entails  the  performance  of  three  (3)  or  more  self-consistent 
calculations.  The  first  calculation  employs  the  minimal  basis  set  for  the  system  and  progressively 
augment  it  in  subsequent  calculations.  The  occupied  energy  levels  or  bands  of  a  calculation  the 
next,  consecutive  calculation  whose  basis  set  is  larger.  As  explained  in  detail  in  the  articles, 
successive  calculations  have  different  basis  sets:  the  basis  set  for  the  latter  is  obtained  by 
augmenting  that  of  the  former  with  the  next  excited  state  orbital  an  atomic  or  ionic  species  that  is 
present  in  the  system  under  study  (i.e.,  In  and  N  for  InN).  This  process  continues  until  the 
occupied  energies  are  found  to  be  identical  for  two  consecutive  calculation  N  and  (N+1).  Then, 
the  results  of  calculation  N  are  the  physical  ones  for  the  system  under  study  (and  not  those  of 
calculation  N+1).  This  criterion  for  stopping  the  process  rests  trivially  on  the  facts  that  (a) 
density  functional  theory  is  a  ground  state  theory  and  (b)  only  the  wave  functions  of  the  occupied 
states  are  including  in  the  construction  and  reconstruction  of  the  charge  density  (and  hence  the 
potential  and  the  Hamiltonian)  in  the  self-consistency  process. 

Once  the  correct  energy  levels  or  bands  are  obtained,  then  we  proceed  to  calculate  a  host 
of  properties  that  include  the  energy  or  band  gaps,  optical  transition  energies,  the  total  energy, 
the  total  and  partial  densities  of  states,  the  effective  masses  of  charge  carriers  (for  application  in 
transport  studies),  and  the  optical  properties  (i.e.,  the  complex  dielectric  functions).  Of  course, 
the  total  energy  curves  are  utilized  to  compute  the  bulk  modulus. 

PROJECT  MANAGEMENT  AND  OVERSIGHT  ;  In  addition  to  the  above  technical 
approach,  we  should  note  that  our  programmatic  implementation  also  adhered  to  the  proposal.  In 
particular,  we  have  hired  Dr.  Hua  Jin  as  the  Postdoctoral  Associate.  She  has  been  performing 
excellently,  not  only  in  carrying  out  research  tasks,  but  also  in  assisting  the  investigators  in  their 
training  of  students  on  the  complex  software  package  utilized  in  this  work.  A  graduate  and 
undergraduate  students,  mostly  African  American,  have  also  been  employed.  SUBR  is  a 
Historically  Black  College  and  University  (HBCU).  The  effort  levels  of  the  investigators  (D. 
Bagayoko  and  G.  L.  Zhao)  are  as  per  the  proposal. 

We  timely  submitted  a  progress  report  in  2005,  following  the  instructions  in  the  Award 
Letter,  These  instructions  included  an  address  for  the  Technical  Monitor,  Dr.  Colin  Wood. 
While  the  various  copies  of  the  report  were  received  by  the  concerned  ONR  officials,  the  one  to 
the  Technical  Monitor  was  not  delivered  and  we  did  not  know  this  fact  until  around  the  end  of 
the  summer  of  2006.  [The  address  in  the  Award  Letter  is  totally  different  from  the  one  Dr.  Wood 
provided  to  us  at  the  end  of  the  summer  of  2005.]  We  then  sent  the  complete  report  to  Dr.  Colin 
Wood  as  an  email  attachment  (PDF).  The  second  year  funding  we  were  supposed  to  receive  in 
the  fall  of  2005  was  not  received  until  the  end  of  the  summer  of  2006,  a  month  or  so  from  the 
initial  expiration  date  of  the  project  (in  September  2006).  We  wrote  to  request  a  no-cost 
extension  for  one  year.  At  the  recommendation  of  Dr.  Wood,  we  received  six  month  extension  so 
that  the  actual  expiration  date  of  the  project  is  December  31,  2006.  This  situation  led  to  serious 
strains  in  the  final  phase  of  the  project  where  much  work  had  to  be  done  within  approximately 
four  (4)  months  as  opposed  to  close  to  a  year.  We  managed,  in  addition  to  Dr.  H.  Jin,  to  hire  Ms. 
L.  Franklin  (who  co-authored  2  of  the  attached  refereed  articles),  and  Dr.  Yang  for  a  few  months. 
In  so  doing,  we  completed  the  work.  We  could  not,  however,  hire  a  graduate  student  for  half  a 
year  (with  no  support  in  the  spring  of  2007).  Hence,  the  $12,000  or  so  of  graduate  assistantship. 
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with  the  associated  overhead,  will  be  returned  to  ONR.  The  circumstances  above,  we  hope,  show 
that  these  residuals  were  not  due  to  our  failure  to  manage  the  project  very  well. 


ACCOMPLISHMENTS ;  As  summarized  in  the  abstract  of  this  report,  we  have  published 
twelve  (12)  articles,  made  fourteen  (14)  technical  presentations,  three  of  which  were  at 
international  conferences.  To  reduce  redundancy,  we  refer  the  reader  to  the  following  full  listing 
of  the  articles,  with  referencing  information.  We  also  list  the  presentations  below.  The  subtitles 
below,  for  accomplishments,  are  keyed  to  the  objectives  enunciated  above. 


Publications  (12) 

"Structural,  Elastic,  and  Electronic  Properties  of  Deformed  Carbon  Nanotubes  under  Uniaxial 
Strain.  ”A.  Pullen,  G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang.  Physical  Review  B  71,  205410. 

“Ab-initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short  Carbon  Nanobells.  ”  G. 
L.  Zhao,  D.  Bagayoko,  and  E.  G.  Wang.  Accepted  for  publication  in  the  Proceedings  of  the  2005 
China  Conference  on  Nanoscience  and  Technology,  Beijing,  China,  July  2005. 

“A  Universal  Relation  Between  the  Densities  of  States  Near  van  Hove  Singularities  and  the 
Effective  Electron  Masses  in  1 -Dimensional  Semiconductors.  ”  G.  L.  Zhao  and  D.  Bagayoko. 
Accepted  for  publication  in  the  Proceedings  of  the  2005  China  Conference  on  Nano-Science  and 
Technology  ,  Beijing,  China,  July  2005. 

"LDA  and  LCAO-BZW  Description  of  Electronic  Properties  ofWurtzite  Zinc  Oxide  (w-ZnO).  ” 
Diola  Bagayoko,  Lashounda  Franklin,  and  G.  L.  Zhao,  accepted  for  publications  in  Proceedings 
of  the  2005  National  Conference  of  the  National  Society  of  Black  Physicists  ('www.nsbp.org). 
Orlando,  Florida. 

"Predictions  of  Electronic,  Structural,  and  Elastic  Properties  of  Cubic  InN.  ”  Diola  Bagayoko, 
Lashounda  Franklin,  and  G.  L.  Zhao,  Journal  of  Applied  Physics  96,  4297-4301,  2004. 

‘‘Density  Functional  Band  Gap  of  Wurtzite  InN  ”  Diola  Bagayoko  and  Lashounda  Franklin, 
Journal  of  Applied  Physics,  97,  123708,  2005. 

“Re-examination  of  the  Ab-initio  Calculation  of  the  Electronic  Structures  of  ZnSe,  Ge,  and 
GoAs.  ”  G.  L.  Zhao,  L.  Franklin,  and  D.  Bagayoko.  Submitted  to  Physical  Review  B,  2005. 


Presentations  (14) 

December  13,  2006.  Department  of  Physics,  Cairo  University,  Cairo,  Egypt.  “Predictive 
Calculations  of  Electronic  and  Related  Properties  of  Materials:  The  BZW  Method.  D.  Bagayoko 
(presenter).  Audience:  Condensed  matter  researchers  and  graduate  students. 


June  2,  2006,  Department  of  Physics,  Beijing  (Peking)  University,  Beijing,  China.  “Ab-initio 
Predictions  of  electronic  Properties  of  Materials:  The  Bagayoko,  Zhao,  and  Williams  (BZW) 
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Method.  ”  D.  Bagayoko  (presenter)  and  G.  L.  Zhao.  Audience:  Researchers  and  graduate 
students,  Department  of  Physics,  University  of  Beijing. 

June  9-11,  2005,  Beijing,  China.  China  International  Conference  on  Nanoscience  and 
Technology.  “Ab-initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short  Carbon 
Nanobells.  ”  G.  L.  Zhao,  D.  Bagayoko,  and  E.  G.  Wang. 

June  9-11,  2005,  Beijing,  China.  China  International  Conference  on  Nanoscience  and 
Technology.  “A  Universal  Relation  Between  the  Densities  of  States  Near  van  Hove  Singularities 
and  the  Effective  Electron  Masses  in  1 -Dimensional  Semiconductors.  ”  G.  L.  Zhao  and  D. 
Bagayoko. 

March  15,  2006,  APS  March  Meeting,  Baltimore,  Maryland,  USA.  “Density  Functional  Band 
Gap  ofAlAs.  ”  Hua  Jin,  Guanglin  Zhao,  and  Diola  Bagayoko. 

March  15,  2006,  APS  March  Meeting,  Baltimore,  Maryland,  USA.  “Calculated,  Optical 
Properties  ofWurtzite  InN. "  D.  Bagayoko,  H.  Jin,  and  G.  L.  Zhao. 

March  15,  2006,  APS  March  Meeting,  Baltimore,  Maryland,  USA.  “Predictions  of  Electronic 
and  Optical  Properties  of  Cubic  InN.  ”  L.  Franklin,  H.  Jin,  G.  L.Zhao,  and  D.  Bagayoko. 

March  15,  2006,  APS  March  Meeting,  Baltimore,  Maryland,  USA.  “Optical  Properties  of 
Aligned  Carbon  Nanotube  Mats  for  Photonic  Applications.  ”  G.  L.  Zhao,  D.  Bagayoko,  and  L. 
Yang. 

March  14,  2006,  APS  March  Meeting,  Baltimore,  Maryland,  USA.  “Ab-initio  Simulations  of  the 
Growth  of  Short  Carbon  Nanobells.  ”  G.  L.  Zhao,  and  D.  Bagayoko. 

March  25,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“Structural,  Elastic,  and  Electronic  Properties  of  Deformed  Carbon  Nanotubes  under  Uniaxial 
Strain.  ”  A.  Pullen  SUBR  &  (Caltech),  G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang  (NASA),  Bull. 
APS,  Vol.  50,  No.  1,  Page  1420  (2005). 

March  23,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“Re-examination  of  Ab-initio  Calculation  of  the  Electronic  Structure  of  Zn  Se,  Ge,  and  GoAs. " 
G.  L.  Zhao,  L.  Franklin,  and  D.  Bagayoko,  Bull.  APS,  Vol.  50,  No.  1,  Page  1073  (2005). 

March  22,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“True  EDA  Band  Gaps  of  Wurtzite  and  Cubic  Indium  Nitride  (w-InN  and  c-InN).  "  D.  Bagayoko, 
G.  L.  Zhao,  and  L.  Franklin.  Bull.  APS,  Vol.  50,  No.  1,  Page  617  (2005).  Audience: 
Approximately  70  physical  science  researchers  and  graduate  students. 

February  18,  2005.  Orlando,  Florida,  Disneyland.  2005  National  Conference  of  the  National 
Society  of  Black  Physicists  (NSBP)  and  of  the  National  Society  of  Hispanic  Physicists  (NSHPj. 
“Local  Density  Functional  Description  of  Electronic  Properties  of  Wurtzite  Zinc  Oxide  (ZnO).  ” 
D.  Bagayoko,  G.  L.  Zhao,  and  L.  Franklin.  Audience:  17  faculty  members,  graduate  students, 
and  federal  lab  researchers. 
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February  17,  2005.  Orlando,  Florida,  Disneyland.  2005  National  Conference  of  the  National 
Society  of  Black  Physicists  (NSBP)  and  of  the  National  Society  of  Hispanic  Physicists  (NSHP). 
“A  Competitive  Edge  for  Recruitment:  The  Versatility  and  Wonders  of  Physics.  ”  D.  Bagayoko. 
Audience:  20  faculty  members,  graduate  students,  and  federal  lab  researchers. 

Code  Development 

Another  major  component  of  our  accomplishments  in  this  reporting  period  stems  from  the 
development  and  refinement  of  the  computer  program  package  for  wider  utilization,  including 
accurate  total  energy  calculations.  We  have  also  developed  a  precise  and  concise  user-manual 
that  we  are  currently  field-testing  for  effectiveness. 

Research  Training  of  Students 

We  have  trained  four  (4)  minority  undergraduate  students  in  the  performance  of  the  sophisticated 
calculations  this  project  entails.  These  students  are  physics  majors  at  SUBR.  One  of  the  students 
is  the  first  author  on  the  Physical  Review  article  on  electronic,  elastic,  and  related  properties  of 
single  walled  carbon  nanotubes. 

CONCLUSIONS:  We  have  implemented  in  a  successful  manner,  as  gauged  by  the  publications 
and  technical  presentations,  this  project  focused  on  “predictive  calculations  of  properties  of 
wide-gap  and  nano-semiconductors.  Materials  studied  to  date  include  GaN,  ZnO,  wurtzite  InN, 
cubic  InN,  and  several  carbon  nanotubes.  Our  results,  in  contrast  to  previous  DFT  or  LDA 
calculations  by  other  groups,  are  in  excellent  agreement  with  experiment.  The  significance  of 
these  results  is  underscored  below  and  in  the  abstract. 

SIGNIFICANCE:  Our  articles  on  electronic,  structural,  optical,  and  other  properties  of 
semiconductors  and  carbon  nanotubes  (Please  see  the  attached  list  of  publications)  have 
profound  implications,  even  for  the  study  of  nuclei.  These  papers  resolved  a  problem  dating  back 
to  the  beginning  of  quantum.  Indeed,  we  solved  the  woeful,  theoretical  underestimation  of  the 
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method  of  solution  also  applies  to  the  prediction  of  nuclear  energy  levels  in  the  shell  model, 
opening  the  way  to  theoretical  explorations  of  possibilities  of  a  population  inversion  and  the 
actual  construction  of  a  “gamma  ray  amplification  by  stimulated  emission  of  radiation  (graser)  ” 
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Ab-Initio  Density  Functional  Calculations  of  the  Growth  and 
Structural  Properties  of  Short  Carbon  Nanobells 
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Abstract.  We  performed  ab-initio  density  functional  calculations  to  study  the  structural  and 
growth  properties  of  short  carbon  nanobells.  We  used  a  real  space  approach  and  the  linear 
combination  of  atomic  orbitals  (LCAO)  formalism.  In  the  nitrogen-doped  carbon  nanobells, 
the  nitrogen  atoms  that  are  attracted  to  the  open-edge  sites  of  the  carbon  nanobells  play  an 
important  role  in  the  growth  of  the  short  carbon  nanostructures.  We  present  the  calculated 
electronic  structure  of  the  short  nanobells.  The  calculated  local  density  of  states  of  the 
nanobells  revealed  field  emission  characteristics  that  agree  with  experimental  observations. 


1.  Introduction 

Carbon  nanostructures  present  interesting  properties  with  a  great  potential  of  applications. 
Particularly,  short  carbon  nanotubes  or  carbon  nanobells  present  excellent  field-emission 
properties  that  have  attracted  a  great  interest  for  applications.^’’^^  Experiments  demonstrated 
that  the  nitrogen-assisted  synthesis  can  grow  carbon  nanofibers  on  a  large  scale.^^  '^^  Such 
carbon  nanofibers  exhibit  a  "bamboo-like"  structure.  Distinctively,  a  great  part  of  the 
bamboo-like  nanofibers  consists  of  short  carbon  nanobells.  Individual  nanobells  are  self- 
contained  and  stack  one  on  top  of  the  other  to  create  a  long  nanofiber.  The  carbon  nanobells 
may  be  viewed  as  short  nanotubes  such  that  their  lengths  are  of  the  same  order  as  their 
diameters.  The  electronic  structure  of  the  short  carbon  nanobells  is  substantially  different 
from  that  of  the  pure  and  long  carbon  nanotubes  because  of  the  quantum  effects  inherent  to 
their  sizes.  The  nanobells  exhibit  novel  electron  field  emission  properties  with  a  tum-on  field 
of  electron  emission  as  low  as  ! fmi  We  performed  ab-initio  density  functional 
calculations,  aiming  to  understand  the  growth,  structural,  and  electronic  properties  of  short 
carbon  nanobells. 

2.  Method 

Although  the  experimental  synthesis  of  carbon  nanobells  produced  samples  of  various  sizes, 
measurements  revealed  that  the  atomic  structure  of  the  bamboo-like  morphology  in  thick 
nanobells  is  similar  to  that  of  the  thinner  ones.^^^  Experiments  also  revealed  that  adjacent 
nanobells  do  not  have  a  firm  contact,  but  instead  the  closed  end  of  one  nanobell  is  weakly 
inserted  into  the  open  end  of  another.  Single  nanobells  can  be  easily  separated  from  other 


parts  of  the  nanofiberJ'®^  Hence,  our  ab-initio  calculations  will  focus  on  single  nanobells.  The 
weak  interaction  between  adjacent  carbon  nanobells  may  be  included  as  a  perturbation  in 
further  studies  in  modeling  the  bamboo-like  morphology  of  the  nanofibers.  We  constructed  a 
prototypical  carbon  nanobell  as  a  model  to  simulate  the  general  structure  of  short  carbon 
nanobells.  Figure  1  presents  the  calculated  images  of  the  electron  density  distribution  of  a 
prototypical  carbon  nanobell  in  real  space.  The  prototypical  carbon  nanobell  includes  150 
atoms  and  consists  of  half  of  the  C240  fullerene  with  30  extra  atoms  on  the  bell.  Its  length 
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and  diameter  are  8.5  and  7.4  A ,  respectively.  Although  the  size  of  this  prototypical  carbon 
nanobell  may  differ  from  that  of  fabricated  samples,  their  fundamental  features  are  similar. 
These  features  include  the  structural  properties  of  the  open-edge,  the  wall,  and  the  closed  cap 
of  the  nanobells. 


Figure  1.  The  calculated  images  of  the  electron  density  distribution  in  real  space  for  a  prototypical 
carbon  nanobell  that  is  viewed  from  different  directions 

Our  ab-initio  calculations  are  based  on  the  density  functional  theory  of  Hohenberg- 
Kohn  and  Kohn-Sham.^''*  '^  '^  '’’'^^  In  the  linear  combination  of  atomic  orbitals  (LCAO) 
method,  we  solved  the  Kohn-Sham  equations  self-consistently  by  employing  the  Rayleigh- 
Ritz  variational  process.  In  the  self-consistent  calculations  for  the  electronic  structure 

of  short  CNx  nanobells,  we  employed  an  extended  basis  set  that  includes  atomic  orbitals  of 
C(ls2s3s  2p3p)  and  N(ls2s3s  2p3p).  Here  C(3s  3p)  and  N(3s  3p)  are  extra  orbitals  that  are 
used  to  augment  the  basis  set  to  account  for  possible  charge  diffusion  and  polarization  in  the 
short  nanobells.  The  real  space  approach  of  the  LCAO  method  enables  us  to  complete  the 
required  computations  using  our  SiliconGraphics  0rigin2000  that  is  equipped  with  4  GB 
RAM  (memory). 

3.  Results 

3.1 .  Relaxed  Atomic  Structure  near  the  Open-Edge  of  the  Short  Carbon  Nanobell 
The  open-edges  of  carbon  nanobells  play  an  important  role  in  their  electronic  properties.  The 
structure  near  the  open  edge  of  the  carbon  nanobell  is  different  from  that  in  the  inner  wall  and 
the  closed  cap.  We  performed  total  energy  minimization  to  identify  the  atomic  structure  near 
the  open-edge  of  the  prototypical  carbon  nanobell.  The  calculated  C-C  bond  length  on  the 


first  atomic  ring  at  the  open-edge  of  the  carbon  nanobell  is  1 .37  ^ ,  which  is  much  shorter  than 

the  C-C  bond  length  of  1.415  yi  in  the  inner  atomic  rings  of  the  bell.  This  reduction  of  the  C- 
C  bond  length  for  the  carbon  atoms  at  the  open-edge  sites  is  a  result  of  the  reduction  in 
coordination  number.  The  calculated  bond  length  of  the  carbon  atoms  on  the  first  atomic  ring 
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to  the  carbon  atoms  on  the  second  ring  is  1 .40^4 ,  which  is  not  much  different  from  that  in  the 
inner  rings,  away  from  the  open-edge  of  the  nanobell.  This  result  indicates  that  the  effect  of 
further  relaxation  of  the  third  atomic  ring  and  inner  rings  would  be  insignificant. 

3.2.  A  Possible  Growth  Mechanism  of  Carbon  Nanobells 

The  growth  mechanism  of  the  carbon  nanobells  remains  a  difficult  problem.  There  are 
probably  two  major  growth  models.  The  first  one  assumes  that  the  carbon  atoms  are  added 
at  the  open-ends  of  the  nanobells. The  second  one  involves  the  C2  absorption  process  that 
is  assisted  by  the  pentagonal  defects  on  the  closed  caps.  Although  both  models  are  very 
interesting,  we  have  technical  difficulties  using  density  functional  computations  to  simulate 
the  second  growth  model.  We  successfully  performed  local  density  functional  computations 
to  test  the  first  growth  mechanism.  Because  of  the  active  dangling  bonds  of  carbon  atoms  at 
the  open-edge  of  nanobells,  there  is  a  high  possibility  that  carbon  atoms  can  be  attracted  to 
these  sites  for  growth  to  occur.  We  calculated  the  total  energies  in  the  following  two  cases.  In 
the  first  case,  a  ring  of  20  carbon  atoms  attaches  to  the  open-edge  of  the  carbon  nanobell  at  a 
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bond-length  of  1.415  ^4 .  In  the  second  case,  these  20  carbon  atoms  are  free.  The  calculated 
total  energy  in  case  I  is  substantially  lower  than  the  corresponding  value  in  case  II.  The  total 
energy  difference,  which  may  also  be  defined  as  the  cohesive  energy  of  the  carbon  atoms  in 
the  first  growth  model,  is  4.7  eV/atom.  In  this  growth  model,  see  Figure  1,  C2  dimers  are 
absorbed  at  the  active  dangling  bond  sites  at  the  open  end  of  the  carbon  nanobell.  A  C2  dimer 
that  deposits  on  the  open-edge  of  the  bell  forms  one  covalent  bond  on  the  same  atomic  ring. 
Each  of  the  C  atoms  also  forms  one  covalent  bond  with  the  C  atom  of  the  next  atomic  ring 
and  remains  an  active  dangling  bond  towards  the  open-space.  This  growth  model  is  consistent 
with  previous  analysis  of  experimental  results.^”’'^’^’^  Although  we  cannot  rule  out  the  second 
growth  model  (growth  on  the  closed  cap  of  carbon  nanobells),  our  calculations  presented 
clear  evidence  that  the  first  growth  mechanism  is  highly  possible. 

3.3.  Nitrogen-Doped  Carbon  Nanobell 

Experimental  results  indicated  that  the  growth  of  the  short  carbon  nanobells  is  highly 
dependent  on  the  nitrogen  concentration  in  the  gas  mixture  during  the  synthesis.^^’’®^  Without 
nitrogen  atoms  in  the  growth  gas  mixture  under  the  same  conditions,  long  and  pure  carbon 
nanotubes,  without  the  “bamboo-like”  morphology,  are  produced.  However,  it  was  not  clear 
why  the  nitrogen  atoms  could  turn  the  growth  of  the  would-be  long  carbon  nanotubes  into  that 
of  short  carbon  nanobells.  We  performed  total  energy  calculations  to  study  the  nitrogen- 
doped  carbon  nanobells  (CNx  nanobells).  Since  nitrogen  atoms  can  also  form  the  planar  sp^- 
hybrid,^^^^  we  studied  the  substitutional  doping  of  nitrogen  atoms  in  the  CNx  nanobells.  The 
substitutional  doping  of  nitrogen  atoms  in  CNx  nanobells  was  also  proposed  from  the  analysis 
of  experimental  results.^'^^ 

We  first  replaced  carbon  atoms  with  nitrogen  atoms  at  the  open-edge  of  the 
prototypical  nanobell.  We  compared  the  total  energies  in  the  two  cases.  In  the  first  case, 
computations  are  carried  out  for  the  pristine  carbon  nanobell  in  the  presence  of  ten  free  atoms 


of  nitrogen.  In  the  second  case,  ten  nitrogen  atoms  replace  ten  carbon  atoms  on  the  open-edge 
and  the  substitutionally  N-doped  nanobell  is  in  the  presence  of  ten  free  atoms  of  carbon.  The 
calculated  total  energy  in  the  second  case  is  lower  than  that  in  the  first  case  by  an  amount  of 
0.5  eV  per  atom,  without  relaxing  the  structure.  We  then  performed  total  energy  calculations 
and  identified  the  atomic  positions  of  the  relaxed  structure  of  the  nitrogen-doped  nanobell 
from  the  total  energy  minimization.  The  calculated  C-N  bond  length  in  the  first  atomic  ring  at 

the  open-edge  of  the  nanobell  is  1 .376  A .  The  bond-length,  from  the  nitrogen  atoms  of  the 

first  atomic  ring  to  the  carbon  atoms  of  the  second  atomic  ring,  is  1 .377  A  .  They  are  close  to 

the  length  of  the  partially  double  C-N  bond  (1.352^)  in  heterocyclic  systems.^^^^  The  C-C 

bond  length  from  the  first  atomic  ring  to  the  second  ring  is  1 .397  A  .  The  total  energy  of  the 
nitrogen-doped  nanobell,  in  the  relaxed  structure,  is  lower  than  that  of  the  undoped  one  by 
0.72  eV  per  atom. 

We  also  studied  the  substitutional  doping  with  nitrogen  atoms  (for  carbon  atoms)  in 
the  wall  (away  from  the  open-edge  and  the  closed  cap)  of  the  prototypical  carbon  nanobell. 
We  used  the  same  computational  procedure  as  discussed  above.  The  calculated  total  energy  of 
the  N-doped  nanobell,  at  the  wall  position,  is  much  higher  than  that  of  the  undoped  bell  by  an 
amount  of  5  eV  per  atom.  This  result  indicates  that  the  substitutional  doping  with  nitrogen  in 
the  wall  of  carbon  nanobells  is  not  energetically  favored. 

These  calculated  results  for  the  nitrogen-doped  nanobells  in  the  prototypical  model 
thus  indicate  that  the  dopant  nitrogen  atoms  prefer  to  stay  at  the  open  edge  of  the  nanobell  as 
opposed  to  being  in  the  wall  (i.e.  lower  rings)  of  the  nanobell.  Once  there  are  enough 
nitrogen  atoms  on  the  open  edge  of  the  nanobell,  carbon  atoms  cannot  attach  to  the  open-edge 
of  the  CNx  nanobell.  Such  an  attachment  will  result  in  nitrogen  atoms  being  in  inner  rings 
other  than  that  at  the  open-edge,  resulting  in  a  configuration  that  is  not  favored  energetically. 
Consequently,  nitrogen  atoms  act  as  the  stopper  of  the  growth  of  the  carbon  nanobell.  This 
growth  mechanism  for  the  nitrogen-doped  carbon  nanobells  and  the  effect  of  the  nitrogen 
atoms  in  the  growth  process  agree  with  the  observations  from  experiments.^^’*'’’’ 

We  present  the  calculated  local  density  of  states  (LDOS)  of  the  CNx  nanobell  in 
Figure  2.  The  dopant  nitrogen  atoms  stay  at  the  open-edge  sites  in  the  prototypical  nanobell. 
The  LDOS  from  the  nitrogen  and  carbon  atoms  on  the  first  atomic  ring  at  the  open-edge  of 
the  nanobell  are  shown  in  the  first  and  the  second  panels  of  Figure  2,  respectively.  There  is 
no  N  atom  on  other  atomic  rings  away  from  the  open-edge  of  the  nanobell.  Figure  2  clearly 
shows  that  the  N  and  C  atoms  on  the  first  atomic  ring  at  the  open-edge  of  the  nanobell  have 
the  dominant  contribution  to  the  electronic  states  around  the  Fermi  level.  The  LDOS  at  the 
Fermi  level,  from  the  N  and  C  atoms  on  the  first  atomic  ring  at  the  open-edge,  are  0.56  and 
0.71  states  per  eV  per  atom,  respectively.  The  contribution,  from  the  second  atomic  ring  near 
the  open-edge,  to  the  LDOS  at  the  Fermi  level  is  0.23  states  per  eV  per  atom.  The  third 
atomic  ring  and  the  inner  rings,  that  include  the  closed  cap,  have  a  much  smaller  contribution 
to  LDOS  at  the  Fermi  level,  at  about  0.08  states  per  eV  per  atom.  Thus,  electrons  will  be 
emitted  from  the  atomic  sites  near  the  open-edge  of  the  nanobell  in  field-emission 
experiments.  The  real  nanobell  samples  in  experimental  studies  mostly  involve  multilayers  of 
graphite  sheets  and  the  nanobells  stacked  one  on  top  of  the  other.  The  weak  interactions 
between  graphite  layers  in  these  nanobells  and  that  between  the  bells  may  further  broaden  the 
electronic  energy  levels. 
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Figure  2.  The  calculated  local  density  of  states  of  electrons  of  a  nitrogen-doped  carbon  nanobell.  The 
1®'  ring-N  and  -C  refer  to  the  nitrogen  and  carbon  atoms  on  the  first  atomic  ring  at  the  open-edge  of 
the  nanobell.  The  2"'*  and  3'^“'  rings  refer  to  the  second  and  third  atomic  rings  near  the  open-edge.  The 
inner  rings  refer  to  the  atomic  sites  away  from  the  open-edge,  including  the  closed  cap.  The  Fermi 
level  is  at  0.0  eV. 

4.  Conclusion 

In  conclusion,  we  performed  ab-initio  density  functional  calculations  to  simulate  the  growth, 
structural,  and  electronic  properties  of  prototypical  carbon  nanobells.  In  the  nitrogen-doped 
carbon  nanobells,  nitrogen  atoms  that  are  attracted  to  the  open-edge  of  the  nanobells  may  play 
a  role  of  stopper  of  the  growth  of  the  nanostructures.  The  calculated  local  densities  of  states  of 
the  CNx  nanobells  indicate  that  electrons  are  most  likely  emitted  from  the  atomic  sites  near  the 
open-edge  of  the  carbon  nanobell  in  field-emission  experiments.  This  result  agrees  with 
experimental  observations. 
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Abstract 

We  report  ab-initio,  self-consistent  results  for  the  dielectric  function  of  wurtzite  indium  Nitride 
(w-InN).  Our  calculations  employed  a  local  density  approximation  (LDA)  potential,  the  linear 
combination  of  atomic  orbitals  (LCAO),  and  the  Bagayoko,  Zhao,  and  Williams  (BZW)  method. 
Our  findings  agree  very  well  with  recent  measurements  up  to  photon  energies  of  6  eV.  This 
excellent  agreement  shows  the  correct  description,  by  the  LDA-BZW  method,  of  the  relative 
separations  between  upper  valence  bands  and  low-lying  conduction  bands,  in  general,  and 
corroborates  our  previous  result  of  0.88  eV  for  the  intrinsic,  fundamental  band  gap  of  w-InN,  in 
particular.  We  also  report  results  of  simulations  of  the  effect  of  high  electron  doping  on  the  optical 
properties  of  InN. 
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features  below  1  eV  and  exhibit  a  clear  indication  of  an  absorption  edge  around  2  eV,^"^  in 
accordance  with  the  findings  of  the  experimental  Group  I  for  the  band  gap.  In  contrast,  the  second 
set^''°  of  experiments  not  only  shows  the  onset  of  absorption  below  1  eV  but  also  shows  no 
features  around  2  eV.  A  motivation  for  this  work  is  to  resolve  the  above  discrepancy  between  two 
sets  of  experimental  measurements  of  the  dielectric  functions  (i.e.,  imaginary  and  real  parts)  of  the 
w-InN.  The  above  excellent  results  of  the  DFT-BZW  approach,  for  the  band  gap  and  the  effective 
mass,  portend  its  capability  in  resolving  this  second  controversy  surrounding  the  optical  properties 
of  w-InN  by  comparing  the  calculated  dielectric  function,  obtained  with  the  ab-initio, 
self-consistent,  DFT-BZW  bands,  to  experimental  measurements.  While  the  results  of  these 
calculations  of  the  optical  properties  of  w-InN  are  expected  to  show  an  absorption  edge  type 
feature  around  0.8  eV,  it  is  not  a  priori  known  whether  or  not  they  will  also  exhibit  another  feature 
around  2  eV,  as  suggested  by  the  first  set  of  measurements. 

While  the  above  discrepancy  amply  points  to  the  need  for  this  work,  this  need  is  further 
exacerbated  by  the  fact  that  previous  calculations  entailed  either  fittings'^  to  models  of  the 
dielectric  function  or  a  scissors-type  approximation  to  adjust  the  calculated  imaginary  part  of  the 
dielectric  function'^  or  the  calculated  band  gaps.’^'^^  The  modeling  work  of  Djurisic  and  Li'*  and 
the  calculations  of  Abbar  et  al.*^  generally  agreed  with  the  first  set  of  measurements,  i.e.,  with  an 
absorption  onset  around  2  eV,  as  per  the  imaginary  part  of  the  dielectric  function.  Furthmiiller  et 
al.'^  and  Bechstedt  et  al.’^  following  scissors-type  adjustments  to  0.81  eV  of  their  otherwise 
negative,  calculated  band  gap,  reproduced  dielectric  functions  in  general  agreement  with  the 
results  of  the  second  set  of  measurements. 
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in.  Results:  Calculated  Dielectric  Functions 


The  calculated  £2(0))  for  wurtzite  InN  was  resolved  into  two  components,  i.e.,  e2xy(®)  and 
£2z(co)  that  are  the  averages  of  the  spectra  for  polarizations  perpendicular  and  parallel  to  the  c-axis, 
respectively.  Figures  la  and  lb  show  e2xy(®)  and  e2z(co)  ,  respectively,  along  with  recent 
experimental  data  of  Goldhahn  et  al.^  The  peak  positions,  shapes,  and  intensities  of  our  calculated 
£2(0)),  in  both  polarization  directions,  are  found  to  be  in  agreement  with  experiment  up  to  photon 
energies  of  6  eV,  especially  for  the  component  in  the  xy  plane.  The  calculated  dielectric  functions 
exhibit  the  anisotropy  revealed  by  experimental  measurements.  The  first  peaks  at  0.88  eV  ,  for 
£2xy(0)),  and  at  0.85  eV,  for  £2z(0)),  indicate  the  onset  of  band  edge  absorption.  They  are  measures 
of  the  direct  optical  band  gap.  Our  minimum,  calculated  band  gap  is  0.75  eV,  which  is  lower  than 
the  actual  theoretical  minimum  of  0.88  eV^  due  to  constraint  stemming  from  using  the  same 
number  of  orbitals  of  the  respective  s,  p,  and  d  states  of  In  and  N.  Due  to  this  restriction  that  is 
currently  inherent  to  our  optical  property  calculation  program,  the  basis  set  for  N,  in  the 
self-consistent  calculations  performed  here,  is  slightly  larger  than  the  optimal  one.*  This  value  of 
0.75  eV  still  falls  within  the  range  of  0.7  to  1 .0  eV,  as  per  the  experimental  Group  II  noted  above. 
While  the  above  agreement  of  our  calculated  results  with  experiment  is  expected  in  the  vicinity  of 
the  calculated  band  gap,  the  absence  of  features  indicative  of  absorption,  around  2  eV,  is  an 
important  result.  It  partly  settles  the  debate  between  experiments  in  Set  I  and  Set  II  by  establishing 
that  high  quality  w-InN,  in  the  absence  of  high  free  carriers  concentrations,  behave  as  per  the 
findings  of  experiments  in  Set  II.  Calculated  and  experimental  peak  positions  of  the  imaginary  part 
of  the  dielectric  function,  along  with  the  possible  optical  interband  transitions,  are  listed  in  Table  I. 
Besides  the  gap  structure,  £2xy(ff))  shows  six  peaks  which  are  related  to  critical  points  of  the  band 
structure  (CPBS),  while  £2z(co)  exhibits  three  CPBS.  The  possible  optical  transitions 
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shifted  £2(0))  spectra.  A  typical  calculation  result  is  shown  in  Figure  3.  We  depict  the  calculated 
£2(0))  in  the  two  polarization  directions,  along  with  the  results  of  an  experiment^  in  Set  I,  in  panels 
3a  and  3b,  respectively^.  The  optical  absorption  edge  of  1.9  eV,  the  previous  band  gap,  was 
obtained  when  we  increased  the  Fermi  level  from  the  top  of  valence  band  (TOV)  to  1.85  eV  above 
the  TOV.  The  solid  lines,  in  Figures  3a  and  3b,  represent  the  results  of  the  simulation.  Experiment 
and  simulation  appear  to  agree  reasonably  for  photon  energies  up  to  5  eV. 

V.  Summary 

In  summary,  our  ab-initio,  self-consistent,  DFT-BZW  calculations  have  practically 
reproduced  the  experimental  values  of  the  dielectric  function  of  wurtzite  InN  for  photon  energies 
up  to  6  eV.  DFT-BZW  calculations  are  fundamentally  for  the  description  of  the  ground  state. 
Hence,  the  disagreement  between  DFT-BZW  results  and  experiments,  for  high  photon  energies, 
was  expected.  Given  that  the  work  functions  of  many  semiconductors  is  around  or  below  6  eV,  the 
above  disagreement  for  high  photon  energies  may  not  be  a  serious  handicap  for  InN.  For  other 
semiconductors,  a  key  issue  consists  of  the  determination  of  the  photon  energy  up  to  which 
DFT-BZW  and  experiment  agree.  For  photon  energy  above  6  eV,  the  CPBS  still  agree  fairly  with 
experiment,  as  per  the  values  in  Table  I.  The  optical  absorption  edge  of  1 .9  eV  was  simulated  by 
increasing  the  Fermi  Level  into  the  conduction  band.  The  agreement  between  the  results  of  this 
simulation  and  optical  properties  measurements  of  w-InN  samples  with  a  large  band  gap  (i.e., 
around  2  eV)  appears  to  confirm  that  the  increased  Fermi  level,  as  a  result  of  high  electron  doping, 
is  responsible  for  the  large  band  gaps. 
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Table  1.  Peak  positions  (in  eV)  of  the  imaginary  part  of  the  dielectric  function  82(a))  for  InN. 


Exp.^ 

Exp." 

Cal.*’ 

Cal.*’ 

Optical 

Energy 

^2xy 

£22 

^2xy 

£22 

Transition 

(eV) 

4.88* 

4.72 

E^.sv  -  r2c 

4.62 

5.35 

5.38 

5.54 

5.54 

Miv  “Mic 

5.64 

6.05* 

6.16 

6.62 

Li,2v-Li,2c 

6.09 

L3,4v'Li_2c 

6.18 

M2v“M]c 

6.52 

■A.5,6v-Ai_2c 

6.57 

7.87 

7.63 

8.33* 

M2V-M2C 

7.94 

Msv— M2C 

8.49 

8.60 

8.85 

8.95 

r4,5v-  Esc 

8.60 

Ai_4v-A3_4c 

8.72 

K,v-K,e 

8.87 

K,v-K2,3. 

8.91 

9.44 

9.51* 

9.90 

L3,4v‘L3_4c 

9.34 

M4V — 

9.52 

Shoulders. 


^Experiment:  Spectroscopic  ellipsometry  on  MBE  grown  InN  film.^ 
'’Our  calculation  results. 
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Abstract 

Nous  rapportons  des  resultats  inedits,  sur  les  bandes  d’energies  de  certains  raateriaux,  obtenus 
avec  la  methode  de  Bagayoko,  Zhao,  et  Williams  (BZW).  Les  actes  du  premier  Symposium 
Malien  des  Sciences  Appliquees,  en  2000,  contierment  une  description  detaillee  de  la  methode 
BZW,  en  firanfais.  Des  references  citees  ci-dessous  decrivent  la  methode  en  anglais.  La  methode 
a  resolu  un  probleme  qui  date  du  debut  du  20ieme  siecle.  Une  motivation  pour  cet  article 
consiste  a  encourager  des  collegues  africains  a  utiliser  la  methode  pour  se  placer  a  la  frontiere  de 
la  description  quantique  des  proprietes  electroniques  des  atomes,  molecules,  et  solides,  avec  une 
capacite  de  predire  les  dites  proprietes,  a  commencer  par  les  fosses  d’energie  ou  de  bandes. 

We  report  new  results,  on  energies  bands  of  some  materials,  obtained  with  the  Bagayoko,  Zhao, 
and  Williams  (BZW)  method.  The  Proceedings  of  the  first  Malian  Symposium  of  Applied 
Sciences  (MSAS),  in  2000,  contain  a  detailed  description  of  the  BZW  method  in  French. 
References  cited  below  (1-11)  describe  the  method  in  English.  The  method  has  resolved  a 
problem  dating  fi-om  the  beginning  of  the  20th  Century.  One  motivation  for  this  article  consists 
of  encouraging  some  African  colleagues  to  utilize  the  method  in  order  to  place  themselves  at  the 
forefi-ont  of  the  quantum  mechanical  description  of  electronic  properties  of  atoms,  molecules, 
and  solids,  with  a  capacity  of  predicting  these  properties,  beginning  with  energy  gaps  and  band 
gaps. 

Introduction 

Les  Actes  du  Symposium  Malien  des  Sciences  Appliquees  de  2000  dorment  une  description 
detaillee  de  la  methode  de  Bagayoko,  Zhao,  et  Williams.  Le  probleme  resolu  par  la  methode 
consistait,  durant  la  plupart  du  20ieme  siecle,  d’un  desaccord  flagrant  entre  les  resultats 
experimentaux  et  les  valeurs  calculees  de  la  hauteur  ou  de  I’epaisseur  des  bandes  interdites  des 
semi-conducteurs.  Cette  hauteur,  pour  un  systeme  donne,  consiste  en  la  difference  entre  le 
niveau  vide  le  plus  bas  et  le  niveau  occupe  le  plus  eleve ;  nous  I’appelons  le  fosse  d’energie 
(pour  les  atomes,  molecules,  et  autres  systemes  finis)  ou  fosse  de  bande  pour  les  semi- 
conducteurs  et  les  isolants.  Ce  probleme,  pendant  un  siecle,  a  empeche  la  theorie  (i.e.,  les  calculs 
quantiques)  de  bien  informer  et  de  guider  la  conception  et  la  fabrication  d’ instruments  et 
d’appareils  electroniques  de  tout  genre.  En  effet,  pour  toute  application  d’un  atome,  d’une 
molecule,  ou  d’un  semi-conducteur,  le  fosse  d’energie  ou  de  bande  est  tres  souvent 
incontoumable.  Par  exemple,  I’absorption  d’une  radiation  electromagnetique  par  un  semi- 
conducteur  demande  que  I’energie  de  la  radiation  soit  egale  ou  plus  grande  que  le  fosse  de 
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Des  Nouveaux  Resultats  de  la  Methode  BZW 


Les  nouveaux  resultats  de  la  methode  BZW  comprennent  les  solutions  des  problemes  illustres 
dans  les  tables  ci-dessus.  En  effet,  dans  les  sources  indiquees  pour  les  donnees  dans  ces  tables, 
Bagayoko  et  collaborateurs  ont  rapporte  des  fosses  d’energies  pour  I’oxide  de  zinc  (ZnO)  et  pour 
InN  qui  sont  en  accord  avec  Texperience.  Dans  le  cas  de  InN,  notre  resultat  de  0.88  eV  est  en 
accord  avec  le  group  2  d’experience.  Nous  avons  trouve  une  valeur  de  3.28  eV  pour  ZnO,  bien 
en  accord  avec  T  experience  et  plus  du  double  de  la  plupart  des  autres  valeurs  theoriques 
ob tenues  sans  la  methode  BZW.  La  table  suivante  montre  d ’autres  resultats  de  la  methode  BZW 
pour  InN,  dans  les  reseaux  wurtzite  et  cubique. 


Band  Gaps  (Fosses  de  Bande) 

Calculated 

Experiment 

Wurtzite  (w-InN),  LDA-BZW 
a=  3,544  A,  c=  5,718  A,  u  = 

0,3790 

0,88  eV 

0,7-l,OeV 

As  per  Exp. 

Group  2 

Wurtzite  (w-InN),  GGA-BZW 
a  =  3,544  A,  c  =  5,718  A,  u  =  0,3790 

0,81  eV 

Zinc  Blende  (c-InN),  LDA-BZW 
a  =5,017  A  (Equilibrium) 

0,65  eV 

Not  Available 

Zinc  Blende  (c-InN),  LDA-BZW 
a  =4,98  A  (Experimental) 

0,74  eV 

Not  Available 

Effective  Mass  (Masse  Effective) 

Directions 

Experiment 

Wurtzite  (w-InN),  LDA-BZW 
a=  3,544  A,  c=  5,718  A,  u  = 

0,3790 

r-M  0,088  mo 
r-A  0,082  mo 
r-  K  0,088  mo 

0,07  estimated 
with  a  band  gap 
of  0,7  eV 

Wurtzite  (w-InN),  GGA-BZW 
a=  3,544  A,  c=  5,718  A,  u  = 

0,3790 

r-M  0,081  mo 
"  -  A  0,077  mo 
r-  K  0,082  mo 

0,07  estimated 
with  a  band  gap 
of  0,7  eV 

Zinc  Blende  (c-InN),  LDA-BZW 
a  =  5,017  A  (Equilibrium)  I 

”  -  L  0,065  mo 

-  X  0,066  mo 

K  0,066  mo 

Slot  Available 

Zinc  Blende  (c-InN),  LDA-BZW  [ 

a  =  4,98  A  (Experimental)  [ 

r 

"-L  0,076  mo  1 

^  -  X  0,073  mo 

K  0,073  mo 

slot  Available 
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s’averent  necessaires  dans  de  nombreux  projets  de  recherche.  Inutile  d’ajouter  que  les  resultats 
de  ces  calculs  peuvent  etre  publies  dans  les  joumaux  les  plus  reputes. 

Nous  estimons  que  la  formation  pour  I’execution  de  ces  calculs  ne  prend  que  quelques  mois. 

Un  seul  chercheur  forme  en  la  matiere  peut  encadrer  plusieurs  etudiantes  (ou  etudiants)  de 
doctorat  en  meme  temps!  Ajoutons  que  les  diverses  applications  soulignees  ci-dessous  ont  trait 
non  seulement  aux  secteurs  economiques,  industriel,  et  autres  (i.e.,  pharmaceutiques),  mais  aussi 
a  ceux  de  la  defense  et  de  la  securite  ou  surete  nationales. 

Finalement,  nous  repetons  la  conclusion  dans  les  Actes  de  MSAS  2000  :  Apres  avoir  resolu 
le  probleme  precite  (la  catastrophe  des  fosses  d’energie  ou  de  bande),  la  methode  BZW  a 
ouvert  la  voie  pour 

•  La  prediction  de  nouvelles  molecules  (larges  ou  petites,  inorganiques  ou  organiques) ; 

•  La  prediction  des  fosses  d'energie  ou  de  bande  (band  gaps)  des  materiaux  (pour 
guider  la  conception  et  la  fabrication  d’instruments  electroniques  divers— des  plus  simples 
aux  plus  complexes) ; 

•  Une  comprehension  approfondie  de  la  structure  des  noyaux  atomiques ; 

•  L’amelioration  des  simulations  en  foumissant  des  parametres  pour  les  formes  d’energie 
potentielle  utilisees  par  les  simulations  de  la  dynamique  moleculaire  (DM  or  MD)  ou  de 
la  methode  de  Monte  Carlo  (MC) ;  et 

•  Un  avancement  fundamental  de  notre  connaissance  de  la  matiere  (nucleaire, 
atomique,  moleculaire,  et  condensee)  et  des  limitations  reelles  de  la  theorie  de  la 
fonctionnelle  de  la  densite  ainsi  que  de  celles  d’autres  theories  qui  sont  supposees  la 
corriger  ou  la  remplacer.  Avant  notre  travail  en  1998,  la  theorie  de  la  fonctionnelle  de  la 
densitd  etait  blSmee  par  tous  pour  le  flagrant  disaccord  entre  les  valeurs  calculees  et  les 
valeurs  mesurees  des  fosses  d’energie  ou  de  bande. 

Priere  de  bien  reflechir  sur  les  implications  immediates  et  futures  de  la  methode  BZW  pour  les 
industries  (des  semi-conducteurs,  du  noyau  atomique,  de  la  biotechnologie,  et  surtout  de  la  nano- 
technologie  naissante),  A  I’echelle  du  nanometre,  les  effets  quantiques  sont  non  seulement 
omnipresents,  mais  aussi  ils  ne  sont  pas  negligeables  ! 
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We  studied  the  optical  properties  of  the  aligned  carbon  nanotube  (16,  0),  (10,  0),  and  (8,  4)  mats  for 
photonic  device  applications.  We  employed  ab  initio  density  functional  potentials  and  utlized  the 
linear  combination  of  atomic  orbital  formalism.  We  calculated  the  electronic  structure  of  the  carbon 
nanotube  mats  and  the  real  and  imaginary  parts  of  the  dielectric  functions  as  functions  of  the  photon 
energy.  The  calculated  dielectric  functions  of  the  aligned  carbon  nanotube  mats  show  a  strong 
anisotropy  when  the  electric  field  of  the  light  is  parallel  or  perpendicular  to  the  tube  axes.  Especially, 
there  are  strong  peaks  in  the  imaginary  part  of  the  dielectric  function  near  the  absorption  edges, 
when  the  electric  field  of  the  light  is  parallel  to  the  carbon  nanotube  axes.  The  unusual  optical 
properties  of  the  semiconducting  carbon  nanotube  mats  present  an  opportunity  for  applications  in 
electro-optical  devices  in  the  infrared  energy  region.  ©  2006  American  Institute  of  Physics. 

[DOI:  10.1063/1.2201738] 


I.  INTRODUCTION 

Carbon  nanotubes  possess  unique  electronic  properties 
that  are  very  useful  for  building  electro-optical  devices  at 
nanometer  scales.  Recently,  Wu  et  al.  reported  ultrathin, 
transparent,  optically  homogeneous,  electrically  conducting 
films  of  pure  single- walled  carbon  nanotubes.*  The  films  ex¬ 
hibit  optical  transmittance  comparable  to  that  of  commercial 
indium  tin  oxide  in  the  visible  spectrum,  but  far  superior 
transmittance  in  the  technologically  relevant  2-5  /tm  infra¬ 
red  spectral  band.  These  characteristics  indicate  broad  appli¬ 
cability  of  the  films  for  electrical  coupling  in  photonic  de¬ 
vices.  Kim  et  al.  reported  highly  polarized  absorption  and 
photoluminescence  of  stretch-aligned  single-walled  carbon 
nanotubes  dispersed  in  gelatin  films.  The  highly  polarized 
absorption  and  photoluminescence  are  attributed  to  interband 
optical  transitions  in  the  single-walled  carbon  nanotubes 
(SWCNTs).  The  realization  of  highly  aligned  and  lumines¬ 
cent  SWCNT  thin  films  should  contribute  to  the  development 
of  SWCNTs  as  optoelectronic  materials.  Li  et  al.  also  re¬ 
ported  polarized  optical  absorption  spectra  of  single-walled 
4  A  carbon  nanotubes  arrayed  in  the  channels  of  an  A1P04-5 
single  crystal.^  The  measured  absorption  spectra  agreed  well 
with  the  ab  initio  calculations  based  on  the  local  density 
functional  approximation.  Carbon  nanotubes  may  also  be 
used  in  optical  limiting,  nonlinear  optical  devices,  and  other 
applications'*"*** 

Guo  et  al.  reported  the  linear  and  nonlinear  optical  prop¬ 
erties  of  individual  carbon  nanotubes,  utilizing  full-potential 
projected  augmented  wave  (PAW)  method.**  The  calcula¬ 
tions  were  based  on  ab  initio  density  functional  calculation 
within  the  local  density  functional  approximation  (LDA).  A 
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supercell  geometry  was  adopted  so  that  the  nanotubes  are 
aligned  in  a  square  array  with  the  closest  distance  between 
adjacent  nanotubes  being  at  least  6  A.  They  performed  test 
calculations  with  larger  intertube  distances  and  no  discem- 
able  differences  were  found.  In  their  work,  they  aimed  to 
study  the  optical  properties  of  isolated  or  nearly  isolated  car¬ 
bon  nanotubes. 

Machon  et  al.  performed  ab  initio  density  functional  cal¬ 
culations  for  the  optical  properties  of  4  A  diameter  single- 
walled  carbon  nanotubes.*^  The  calculated  optical  properties 
confirmed  the  experimental  results  of  Li  et  al.^  for  the  strong 
anisotropy  of  the  optical  response  of  carbon  nanotubes. 

In  optoelectronic  applications,  a  large  number  of  carbon 
nanotubes  will  be  assembled  in  a  desired  form.  As  identified 
in  the  experimental  studies,  single-walled  carbon  nanotubes 
are  self-assembled  in  a  triangular  lattice  in  bundles,  strands, 
or  mats.*^’*'*  The  photonic  devices  based  on  carbon  nano¬ 
tubes  can  take  advantage  of  the  strong  anisotropy  of  optical 
properties  of  SWCNTs  under  polarized  light.  This  can  be 
achieved  by  aligning  the  SWCNTs  in  a  certain  direction  to 
form  carbon  nanotube  mats.  In  this  work,  we  aimed  to  study 
the  optical  properties  of  aligned  carbon  nanotube  mats 
(CNTMs),  using  ab  initio  density  functional  calculations.  We 
utilized  the  computational  method  of  the  linear  combination 
of  atomic  orbital  (LCAO)  formalism.  We  calculated  the  op¬ 
tical  properties  of  various  CNTMs  that  were  constructed 
from  semiconducting  SWCNTs.  As  a  further  pursuit  of  the 
applications  of  nanodevices,  we  demonstrate  in  this  work 
that  the  aligned  carbon  nanotube  mats  will  have  relevant 
properties  for  photonic  applications. 

II.  METHOD 

We  performed  ab  initio  quantum  calculations  that  are 
based  on  the  density  functional  theory  of  Hohenberg-Kohn 
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FIG.  1.  The  electron  energy  bands  of  SWCNT  (16,  0). 


and  Kohn-Sham.’"’”'^  The  many-body  exchange  and  correla¬ 
tion  interactions  of  the  electron  gas  are  described  by  the  den¬ 
sity  functional  potential  in  a  local  density  approximation 
(LDA).  We  used  the  Cepcrley-Alder  form  of  the  exchange- 
correlation  potential.  As  we  reported  in  the  previous  publi¬ 
cations  on  the  calculations  of  the  electronic  structure  of 
SWCNTs,'®  we  also  performed  ah  initio  calculations  using 
the  generalized  gradient  approximation  (GGA)  potential  of 
Perdew  and  Yue.^*  We  compared  the  electron  energy  bands 
of  the  SWCNTs  from  the  LDA  and  GGA  calculations.  The 
average  difference  between  the  energy  bands,  resulting  from 
the  two  potentials,  was  about  1  mRy.  We  utilized  the  LCAO 
method  in  solving  the  Kohn-Sham  equations.  ’  We  em¬ 
ployed  an  extended  basis  set  that  includes  atomic  orbitals  of 
C(1^2^3^  2/73p).  Here  C{3s  3p)  arc  the  extra  orbitals  that 
are  used  to  augment  the  basis  set  to  account  for  possible 
charge  diffusion  and  polarization  in  the  CNTMs.  The  real 
space  approach  of  the  LCAO  method  enables  us  to  complete 
the  required  large-scale  computations  using  our  Silicon- 
Graphics  0rigin2000  that  is  equipped  with  2  GB  random 
access  memory  (RAM). 

As  identified  in  the  experimental  studies,  single-walled 
carbon  nanotubes  are  self-assembled  in  a  triangular  lattice  in 
bundles,  strands,  or  mats.  We  constructed  the  carbon  nano¬ 
tube  mats  by  aligning  the  SWCNTs  in  an  array  of  the  trian¬ 
gular  lattice.  We  vary  the  distances  between  the  SWCNTs  to 
simulate  the  different  samples  of  the  interactions  in  the  mats. 
We  performed  ah  initio  density  functional  calculations  to 
study  the  electronic  structure  of  the  aligned  carbon  nanotube 
mats.  We  also  carried  out  ah  initio  total  energy  calculations 
for  the  carbon  nanotube  mats  at  various  tube-tube  distances. 

For  the  studies  of  the  optical  properties  of  the  aligned 
CNTMs  for  photonic  applications,  we  are  particularly  inter¬ 
ested  in  the  semiconducting  SWCNTs.  We  utilized  the  cal¬ 
culated  electronic  energy  levels  and  related  wave  functions 
to  evaluate  the  dielectric  function  of  the  aligned  CNTMs.  We 
calculated  the  real  and  imaginary  parts  of  the  dielectric  func¬ 
tion  of  CNTMs  as  functions  of  photon  energies.  The  imagi¬ 
nary  part  of  the  dielectric  function  e2(w),  from  the  direct 
interband  transitions,  is  calculated  from  the  Kubo- 
Greenwood  formula,^'* 
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FIG.  2.  The  electron  density  of  states  of  SWCNT  (16,  0). 
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where  hw  is  the  photon  energy;  P  is  the  momentum  operator, 
P=-/ftV;  B^^|  and  'Pij,  are  the  eigenenergy  and  related  eigen- 
wave  function  of  the  /th  state  at  k  point  in  the  Brillouin  zone; 
fi^i  is  the  Fermi  distribution  function;  and  Cl  is  the  volume  of 
the  triangular  lattice  unit  cell.  The  real  part  of  the  dielectric 
function  8](w)  was  found  from  the  well-known  Kramers- 
Kronig  (K-K)  relation. 


III.  RESULTS 

A.  Single-walled  carbon  nanotube  (16,  0)  mat 

We  constructed  the  carbon  nanotube  mat  utilizing 
SWCNT  (16,  0)  as  the  basic  building  block.  The  diameter  of 
individual  SWCNT  (16,  0)  is  1.25  nm,  which  is  comparable 
to  the  ones  that  were  experimentally  studied  by  Kazaoui 
et  alP 

We  calculated  the  electronic  structure  of  individual 
SWCNT  (16,  0).  There  are  64  atoms  in  the  tube  unit  cell  of 
SWCNT  (16,  0).  The  large  number  of  atoms  per  unit  cell  that 
enter  into  the  ab  initio  calculations  presents  some  technical 
challenges.  We  utilized  a  real  space  approach  of  LCAO  cal¬ 
culations.  The  calculated  electron  energy  bands  of  SWCNT 
(16,  0)  are  presented  in  Fig.  1,  where  r= (0,0,0), 
Z=7r/r(0,0, 1),  and  T  is  the  length  of  the  tube  unit  cell. 
Figure  1  shows  a  direct  band  gap  of  0.62  eV  at  the  F  point. 
The  calculated  density  of  states  (DOS)  of  SWCNT  (16,  0)  is 
shown  in  Fig.  2.  The  van  Hove  singularities  at  the  band 
edges  and  the  sharp  structures  in  the  density  of  states  are 
attributed  to  the  one-dimensional  (ID)  tube  structure. 

Subsequently,  we  calculated  the  electronic  structure  of 
the  carbon  nanotube  mat.  In  contrast  to  the  work  of  Guo 
et  a/.,”  we  aligned  the  SWCNTs  (16,  0)  in  a  triangular  lattice 
to  form  the  mat.  Guo  et  al.  used  a  supercell  geometry,  and 
the  nanotubes  were  aligned  in  a  square  array  with  the  closest 
distance  between  the  adjacent  nanotubes  being  at  least  6  A. 
In  this  work,  we  performed  the  ab  initio  total  energy  calcu¬ 
lations  for  the  carbon  nanotube  mats  at  various  tube-tube 
distances.  The  total  energy  minimization  found  the  stable 
tube-tube  distance  at  3.5  A,  which  is  close  to  the  interplanar 
distance  (3.35  A)  of  graphite. 
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FIG,  3.  The  electron  density  of  states  of  carbon  nanotube  (16,  0)  mat  at  the 
tube-tube  distance  of  3,5  A, 

In  Fig.  3,  we  present  the  electron  density  of  states  of  the 
carbon  nanotube  mat,  which  was  constructed  from  SWCNTs 
(16,  0)  and  the  tube-tube  distance  was  at  3.5  A.  The  density 
of  states  of  the  CNTM  in  Fig.  3  retained  the  van  Hove  sin¬ 
gularities.  The  weak  interaction  between  the  nanotubes 
slightly  smeared  the  sharp  peak  structures  in  the  DOS  as 
compared  to  those  of  SWCNTs  in  Fig.  2. 

Utilizing  the  electron  energy  levels  and  related  wave 
functions  from  the  electronic  structure  computations,  we  fur¬ 
ther  calculated  the  optical  properties  of  the  carbon  nanotube 
mat.  We  calculated  the  imaginary  (£2)  part  of  the  dielectric 
function  as  a  function  of  photon  energy.  The  real  part  of  the 
dielectric  function  was  found  from  the  K-K  relation. 
The  calculated  dielectric  function  revealed  a  strong  aniso¬ 
tropy  when  the  electric  field  (£)  of  the  light  is  parallel  or 
perpendicular  to  the  tube  axis.  We  present,  in  Fig.  4,  the 


Energy  (eV) 

FIG.  4.  The  calculated  real  and  imaginary  parts  of  the  dielectric  function, 
e(a))  =  Ei(&))+ie2(w),  of  the  carbon  nanotube  (16,  0)  mat  which  has  the 
tube-tube  distance  of  3.5  A  between  the  SWCNTs  (16,  0). 
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FIG.  5.  The  joint  density  of  states  for  the  aligned  carbon  nanotube  (16,  0) 
mat. 

calculated  dielectric  function  of  the  aligned  nanotube  mat, 
for  the  tube-tube  distance  of  3.5  A  between  the  individual 
SWCNTs  (16,  0).  Figures  4(a)  and  4(b)  present  the  dielectric 
functions  (81,62)  for  the  electric  field  of  the  light  parallel  and 
perpendicular  to  the  tube  axis,  respectively.  When  the  elec¬ 
tric  field  of  the  light  is  parallel  to  the  tube  axis,  the  imaginary 
part  of  the  dielectric  function  82  has  a  strong  peak  at  the 
photon  energy  of  about  0.58  eV,  as  in  Fig.  4(a).  The  height 
of  the  first  peak  of  82  in  Fig.  4(a)  can  reach  a  value  of  about 
60,  which  is  about  a  factor  of  7  higher  than  that  of  the  first 
peak  of  82  (at  0.9  eV)  in  Fig.  4(b).  The  strong  peak  in  82  in 
the  low  photon  energy  region  in  Fig.  4(a)  shows  a  strong 
anisotropy  in  the  optical  properties  of  CNTM  and  is  spec¬ 
tacularly  different  from  that  of  conventional  optical  materi¬ 
als.  The  above  strong  peak  in  82  indicates  that  the  carbon 
nanotube  mats  can  be  used  as  photonic  sensors  in  the  infra¬ 
red  region.  Such  unusual  optical  properties  in  semiconduct¬ 
ing  carbon  nanotube  mats  present  an  opportunity  for  appli¬ 
cations  in  electro-optical  devices. 

We  further  studied  the  physical  factors  for  the  strong 
anisotropy  in  the  optical  properties  of  the  aligned  carbon 
nanotube  mat.  Particularly,  a  strong  anisotropy  is  observed  in 
the  first  peaks  of  the  imaginary  parts  of  the  dielectric  func¬ 
tion  when  the  polarization  of  the  light  is  parallel  or  perpen¬ 
dicular  to  the  tube  axis.  It  is  noted  that  the  first  absorption 
peak  in  82,  in  the  aligned  CNTM,  is  associated  with  the  van 
Hove  singularities  that  are  located  just  above  and  below  the 
Fermi  level.  We  calculated  a  joint  density  of  states  (JDOS) 
that  is  based  on  the  same  calculation  formula  as  that  of  the 
computation  of  82,  but  without  the  inclusion  of  the  electron 
excitation  matrix  elements  in  the  calculation.  The  calculated 
JDOS  of  the  aligned  carbon  nanotube  (16,  0)  mat  is  shown  in 
Fig.  5.  The  first  peak  of  the  JDOS  in  Fig.  5,  at  0.58  eV,  is 
associated  with  the  first  peak  of  82  in  Fig.  4(a)  when  the 
electric  field  of  the  light  is  parallel  to  the  tube  axis.  The 
second  peak  of  the  JDOS  in  Fig.  5,  at  about  0.9  eV,  is  asso¬ 
ciated  with  the  first  peak  of  82  in  Fig.  4(b)  when  the  electric 
field  of  the  light  is  perpendicular  to  the  tube  axis.  The  height 
of  the  first  peak  of  the  JDOS  in  Fig.  5  is  greater  than  that  of 
the  second  peak,  with  a  ratio  of  1.67.  The  electron  excitation 
matrix  elements  that  contribute  to  the  first  peak  of  82  in  Fig. 
4(a)  have  a  much  higher  value  than  those  contributing  to  the 
first  peak  of  82  in  Fig.  4(b).  The  ratio  of  the  excitation  matrix 
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FIG.  6.  The  calculated  real  and  imaginary  parts  of  the  dielectric  function  of 
the  carbon  nanotube  (10,  0)  mat  which  has  the  tube-tube  distance  of  3.34  A 
between  the  individual  SWCNTs  (10,  0). 

elements  accounting  for  these  two  peaks  is  about  4-5. 
Hence,  the  strong  anisotropy  in  £2  *s  mainly  attributed  to  the 
electron  excitation  matrix  elements  when  the  polarization  of 
the  light  is  in  different  directions. 

B.  Single-walled  carbon  nanotube  (10,  0)  mat 

The  electronic  structure  of  individual  SWCNT  (10,  0), 
from  ah  initio  density  functional  calculations,  was  presented 
in  previous  publications.’^’^®  The  diameter  of  individual 
SWCNT  (10,  0)  is  7.83  A.  In  this  work,  we  calculated  the 
optical  properties  of  the  aligned  carbon  nanotube  (10,  0)  mat 
for  electro-optical  device  applications.  The  SWCNTs  (10,  0) 
are  arranged  in  a  triangular  lattice  in  the  aligned  CNTM.  We 
performed  ah  initio  electronic  structure  calculations  for  the 
aligned  carbon  nanotube  (10,  0)  mats.  We  employed  the 
ah  initio  total  energy  minimization  that  found  a  stable  tube- 
tube  distance  at  3.34  A.  In  Fig.  6,  we  present  the  calculated 
real  and  imaginary  parts  of  the  dielectric  function  of  the  car- , 
bon  nanotube  (10,  0)  mat,  for  the  tube-tube  distance  of 
3.34  A.  The  general  feature  of  the  dielectric  function  of  the 
aligned  carbon  nanotube  (10,  0)  mat  is  similar  to  that  of  the 
(16,  0)  mat.  The  calculated  imaginary  parts  of  the  dielectric 
function  in  Fig.  6(a),  for  the  electric  field  of  the  light  parallel 
to  the  tube  axis,  shows  a  strong  peak  near  the  absorption 
edge  at  about  0.75  eV.  When  the  electric  field  of  the  light  is 
perpendicular  to  the  tube  axis,  the  absorption  is  much 
weaker,  as  shown  in  Fig.  6(b).  This  feature  again  presents  a 
strong  anisotropy  in  the  optical  properties  of  the  aligned  car¬ 
bon  nanotube  (10,  0)  mat.  The  height  of  the  first  peak  in  £2 
Fig.  6(a)  for  the  carbon  nanotube  (10,  0)  mat  is  lower  than 
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FIG.  7.  The  calculated  real  and  imaginary  parts  of  the  dielectric  function  of 
the  carbon  nanotube  (8,4)  mat  which  has  the  tube-tube  distance  of  3.34  A 
between  the  chiral  SWCNTs  (8,  4). 

that  of  the  carbon  nanotube  (16,  0)  mat,  which  could  indicate 
the  effect  of  the  nanotube  curvature  on  the  optical  properties. 

The  calculated,  strong  anisotropy  in  the  optical  proper¬ 
ties  of  the  aligned  carbon  nanotube  (10,  0)  mat  as  well  as  the 
locations  of  the  major  absorption  peaks  in  the  imaginary  part 
of  the  dielectric  function  82  agree  with  those  properties  of  the 
isolated  or  nearly  isolated  carbon  nanotubes  that  were  calcu¬ 
lated  by  Guo  et  al.  using  supercells  of  a  square  array.”  In 
this  work,  we  are  particularly  interested  in  the  aligned  carbon 
nanotube  mats  that  are  constructed  from  semiconducting 
SWCNTs  for  optoelectronic  applications.  Guo  et  al.  pre¬ 
sented  the  ealculated  results  of  ei  and  £2  for  the  isolated 
SWCNT  (10,  0)  (Fig.  8  in  their  article)”  when  the  eleetrie 
field  is  parallel  and  perpendicular  to  the  nanotube  axis.  How¬ 
ever,  they  did  not  present  the  magnitude  of  the  first  absorp¬ 
tion  peak  in  £2  when  the  electric  field  is  parallel  to  the  nano¬ 
tube  axis.  The  location  of  the  first  major  peak  in  £2  for  the 
electric  field  in  parallel  to  the  nanotube  axis  was  found  at 
about  0.75  eV  from  our  calculation.  This  finding  is  in  agree¬ 
ment  with  the  result  of  Guo  et  al.  of  about  0.8  eV,  estimated 
from  their  graph. 

C.  Single-walled  carbon  nanotube  (8,  4)  mat 

The  electronic  structure  of  nonsymmorphic  SWCNT  (8, 
4)  has  been  studied  by  ah  initio  density  functional  calcula¬ 
tions  in  previous  work.’®’^®  The  diameter  of  SWCNT  (8,  4)  is 
8.29  A.  We  constructed  a  carbon  nanotube  (8,  4)  mat  by 
arranging  the  SWCNTs  (8,  4)  in  a  triangular  lattice.  The 
tube-tube  distance  in  the  mat  in  the  triangular  lattice  is 
3.34  A,  as  that  in  the  carbon  nanotube  (10,  0)  mat.  We  per- 
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formed  ab  initio  density  functional  computations  to  study  the 
optical  properties  of  the  carbon  nanotube  (8, 4)  mat.  Figure  7 
presents  the  calculated  optical  properties,  i.e.,  the  real  and 
imaginary  parts  of  the  dielectric  functions  as  functions  of  the 
photon  energy.  Figures  7(a)  and  7(b)  present  the  cases  of  the 
electric  field  polarization  of  the  light  parallel  and  perpen¬ 
dicular  to  the  tube  axis,  respectively.  Similar  to  the  carbon 
nanotube  (16,  0)  and  (10,  0)  mats,  there  is  a  strong  peak  in 
the  imaginary  part  (82)  of  the  dielectric  function  near  the 
absorption  edge  at  about  0.89  eV,  when  the  electric  field  of 
the  light  is  parallel  to  the  tube  axis.  The  imaginary  part  (£2) 
of  the  dielectric  function  can  reach  a  value  of  about  49  at  the 
photon  energy  of  0.89  eV.  There  is  a  strong  anisotropy  in  the 
optical  properties  in  the  carbon  nanotube  (8,  4)  mat,  as 
shown  in  Figs.  7(a)  and  7(b). 

IV.  CONCLUSION 

In  this  work,  we  studied  the  optical  properties  of  the 
aligned  carbon  nanotube  (16,  0),  (10,  0),  and  (8,  4)  mats  for 
photonic  device  applications.  We  employed  ab  initio  density 
functional  computations  to  calculate  the  imaginary  and  real 
parts  of  the  dielectric  functions  as  functions  of  the  photon 
energy.  The  calculated  dielectric  functions  of  the  aligned  car¬ 
bon  nanotube  (16,  0),  (10,  0),  and  (8,  4)  mats  present  a  strong 
anisotropy  when  the  electric  field  of  the  light  is  parallel  or 
perpendicular  to  the  tube  axes.  Especially,  there  are  strong 
peaks  in  the  imaginary  part  of  the  dielectric  function  near  the 
absorption  edges,  when  the  electric  field  of  the  light  is  par¬ 
allel  to  the  carbon  nanotube  axes.  The  unusual  optical  prop¬ 
erties  of  the  semiconducting  carbon  nanotube  mats  present 
an  opportunity  for  applications  in  electro-optical  devices  in 
the  infrared  energy  region. 
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We  present  results  of  ab  initio,  self-consistent  calculations  of  electronic  properties  of  AIAs  in  the  zinc-blende 
structure.  Our  nonrelativistic  calculations  employed  the  generalized  gradient  approximation  of  density  func¬ 
tional  potential  and  Bagayoko,  Zhao,  and  Williams  implementation  method  of  the  linear  combination  of  atomic 
orbitals  formalism.  Our  calculated  indirect  band  gaps  at  the  X  and  L  points  are  2.15  and  2.38  eV,  respectively, 
in  good  agreement  with  experimental  values.  The  calculated  direct  gap  at  F  is  25%  smaller  than  the  experi¬ 
mental  one.  We  also  present  calculated  total  and  partial  densities  of  states  and  the  electron  effective  mass  at  the 
bottom  of  the  conduction  band  at  F. 
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1.  INTRODUCTION 

AIAs  has  attracted  much  attention  because  it  is  not  only 
one  of  the  most  important  electronic  and  optoelectronic  ma¬ 
terials,  but  also  a  very  essential  component  in  GaAs-based 
heterostructures,  supperlattices,  and  quantum  wells.  There 
are  numerous  device  applications  of  Al^^Gai.j^As  alloys,  in¬ 
cluding  diode  lasers,  light-emitting  diodes,  photodetectors, 
and  electro-optic  modulators.  AIAs  and  its  heterostructures 
have  been  studied  extensively  for  both  their  scientific  and 
technological  relevance.’"^ 

The  theoretical  prediction  of  band  structures,  fundamental 
energy  gaps,  and  effective  masses  of  semiconductors  and  al¬ 
loys  is  of  great  importance  for  the  fabrication  of  the  hetero¬ 
structures  and  devices.  Over  the  past  2  decades,  there  have 
been  many  theoretical  calculations  of  the  electronic  struc¬ 
tures  of  AlAs.^’^’’“*®  To  our  knowledge,  however,  except  for 
the  fitting  approaches  such  as  the  tight-binding  model,'^’'"*  '^ 
almost  all  of  the  theoretical  calculations  of  the  band  structure 
of  AIAs  led  to  band  gaps  that  deviate  from  the  experimental 
values^’^”*^  to  varying  degrees.  Table  I  lists  the  referenced 
theoretical  and  experimental  band  gaps  of  AIAs.  For  the  the¬ 
oretical  works,  the  table  also  shows  the  applicable  potential 
and  computational  method.  Unlike  other  III-V  compounds, 
AIAs  is  an  indirect  band  gap  semiconductor  with  the  conduc¬ 
tion  band  minimum  close  to  or  at  the  X  point^”'^  *^  while  the 
valence  band  maximum  is  at  the  gamma  (F)  point.  Usually, 
three  band  gaps  are  reported  for  the  AIAs  ^stem:  the  mini¬ 
mum  indirect  band  gap  the  direct  gap  and  the  second 
and  larger  indirect  gap  As  can  be  seen  from  Table  I, 
previous,  ab  initio  local  density  approximation  (LDA)  calcu¬ 
lations  typically  reported  values  of  the  minimum,  indirect 
gap  that  are  36%  to  47%  off  the  measured  one.'*  Although  the 
agreement  with  experiment  has  been  improved  dramatically 
by  using  the  Green  function  and  screened  Coulomb  potential 
approximation  and  the  pseudopotential  method  (GW-PP),’^ 
the  Green  function  and  screened  Coulomb  potential  approxi¬ 
mation  and  the  quasiparticle  approach  (GW-QP),®  and  a 
semi-a()  initio  approach  utilizing  a  minimum  basis  set  of 
orthogonalized  functions  in  a  linear  combination  of  atomic 
orbitals  (OLCAO)®  and  an  additional  atomiclike  potential, 
the  discrepancies  for  the  calculated  band  gaps  E^,  E^,  and  E^ 
have  remained  significant.  These  band  gaps,  for  the  GW-PP*^ 


and  the  GW-QP  are  2.08,  2.75,  and  2.79  eV  for  the  former 
and  2.09,  3.26,  and  2.81  for  the  latter.  While  the  improve¬ 
ment  for  the  minimum,  indirect  gap  is  significant,  the  over¬ 
estimations  of  the  larger,  indirect  gap  are  by  18%  and  28% 
for  the  GW-PP  and  GW-QP  calculations,  respectively. 

The  above  discrepancies  are  key  motivations  for  this  work 
that  employs  a  GGA  potential  and  the  Bagayoko,  Zhao  and 
Williams  (BZW)  method^®’^*  within  the  linear  combination 
of  atomic  orbital  formalism  (LCAO).  The  BZW  method  has 
been  successfully  applied  to  reproduce  or  to  predict  the  band 
gaps  of  numerous  semiconductors,  including  cubic  InN,^^ 
wurtzite  InN,^^  GaN,  Si,  and  C,^'  and  carbon  nanotubes.^'*’^^ 

II,  COMPUTATIONAL  METHOD  AND  DETAILS 

Our  nonrelativistic  calculations  employed  a  nonlocal  den¬ 
sity  functional  potential  from  the  generalized  gradient  ap¬ 
proximation  (GGA).^^"^®  We  utilized  the  formalism  of  the 
linear  combination  of  atomic  orbitals  (LCAO)  in  real  space. 
The  implementation  of  the  BZW  method  in  carrying  out  the 
self-consistent  computations  is  the  major,  distinctive  feature 
of  the  present  work  as  compared  to  previous  works  on  AIAs. 
The  details  of  the  BZW  procedure  are  widely  available  in  the 
literature^* and  are  discussed  further  below.  A  brief  over¬ 
view  of  its  implementation  follows. 

In  the  implementation  of  the  BZW  procedure,  we  started 
the  self-consistent  calculations  with  a  minimal  basis  set,  i.e., 
the  basis  set  just  accounting  for  all  the  electrons  in  the 
atomic  or  ionic  species  present  in  the  system  under  study. 
For  AIAs,  we  chose  the  ions  Al*'*'  and  As*“  for  our  self- 
consistent  calculations,  as  preliminary  studies  indicated 
charge  transfers  closer  to  these  species,  i.e.,  the  self- 
consistent  system  is  approximately  Al**'^‘*^As*’'^‘*".  We  subse¬ 
quently  carried  out  several  other  self-consistent  calculations 
with  larger  and  larger  basis  sets  by  augmenting  with  one  or 
more  ionic  orbitals  that  belong  to  the  next  and  lowest-lying 
energy  levels  in  Al*'*’  or  As*“.  The  occupied  bands  of  a  given 
calculation  are  compared  to  those  of  the  previous  one  until 
they  are  found  to  be  identical  in  numerical  values,  curvature, 
and  branching,  within  computational  uncertainties.  The  re¬ 
sults  reported  here  are  those  of  the  calculation  before  the  last 
one.  The  basis  set  for  this  calculation  (i.e.,  VII  for  AIAs)  is 
referred  to  as  the  optimal  basis  set.  As  explained  further  be- 
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TABLE  1.  Comparison  of  theoretical  and  experimental  band  gaps  of  zinc-blende  AlAs.  The  three  band 
gaps  at  respective  Tisv-Xic,  ruy-Tic  and  Fisv-Lic. 


Potential 

Computational  method 

(eV) 

e[  (eV) 

E^  (eV) 

Experiment  1 

Photoluminescence  (7=12  K) 

2.25" 

Experiment  2 

Excitonic  gap,  photoluminescence  (7=4  K) 

2.23'’ 

3.13*’ 

Excitonic  gap,  photoluminescence  (7=300  K) 

2.15^’ 

3.03*’ 

Experiment  3 

Transport  (7=295  K) 

2.16‘= 

2.98" 

2.36'= 

GGA 

LCAO-BZW  (Present  work)  a=5.66  A 

2.15 

2.35 

2.38 

GGA 

LCAO-BZW  (Present  work)  a=5.6524  A 

2.14 

2.38 

2.39 

LDA 

Projector-augmented-wave  (PAW) 

1.32‘^ 

1.94“' 

2.06'' 

GW 

PAW 

1.57‘^ 

2.7‘‘ 

2.73  '' 

LDA 

Pseudopotential  method  (PP) 

1.20" 

1.77" 

1.89" 

GW 

PP 

2.08" 

2.75" 

2.79" 

LDA 

Ab  initio  PP 

1.44'' 

2.35f 

2.12f 

LDF 

First-principles  full-potential  self-consistent 
lincarized-muffin-tin-orbital  (LMTO) 

1.318 

Ab-initio  PP 

1.3'’ 

2.5'’ 

LDA 

Self-consistent  full-potential 
lincarized-augmented-plane-wave  (FLAPW) 

1.36‘ 

1.95' 

2.07' 

GW 

Quasiparticle  (QP) 

2.09' 

3.26i 

3.03' 

LDF 

Minimal  basis  semi-a/)  initio  orthogonal i zed 
LCAO  method 

2.37'‘ 

2.79'' 

2.81*' 

Tight-binding 
model  (TBM) 

PP 

2.21' 

2.79' 

2.48' 

TBM 

2.142'" 

2.998"’ 

2.313"’ 

TBM 

2.262" 

2.974" 

2.756" 
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low  in  the  discussion  section,  upon  the  convergence  of  the 
charge  density,  potential,  and  occupied  energy  levels,  as  is 
the  case  for  the  calculation  with  the  optimal  basis  set,  com¬ 
putations  with  larger  basis  set  arc  inherently  affected  by  a 
basis  set  and  variational  effect.  This  effect  stems  from  the 
use  of  a  variational  approach  where  the  charge  density  is 
computed  using  only  the  wave  functions  of  occupied  states 
and  consequences  of  the  Rayleigh  theorem. Table  II 
shows  the  basis  sets  employed  in  the  various,  self-consistent 
calculations  we  performed  for  AlAs.  The  results  discussed 
here  are  those  from  Calculation  VII  for  which  the  basis  set  is 
optimal  in  the  sense  specified  above. 

We  utilized  an  expanded  version  of  the  electronic  struc¬ 
ture  calculation  program  package  from  the  Ames  Laboratory 
of  the  Department  of  Energy  (DOE).^^“^'  AlAs,  in  the  zinc- 
blende  structure,  is  a  member  of  the  III-V  family  and  belongs 
to  the  r^-F43/rj  space  group.  As  per  Table  I,  only  Experi¬ 
ment  3,  done  at  a  temperature  of  295  K,  provides^  values  for 
the  two  indirect  and  one  direct  gaps  noted  above.  For  com¬ 
parison  purposes,  we  selected  the  experimental  lattice  param¬ 
eter  of  a=5.66  A  for  a  temperature  of  291  Except  for 


the  gap  values  shown  in  Table  I  for  a  zero-temperature  lattice 
constant^  of  5.6524  A,  the  results  discussed  here  are  for  the 
above  room  temperature  lattice  constant.  The  atomic  wave 

TABLE  II.  The  atomic  orbitals  used  in  calculations  I  to  VIII  for 
AlAs.  Superscript  zeros  indicate  added  orbitals  representing  unoc¬ 
cupied  atomic  states. 


Basis  set  0;  core-state  orbitals  in  calculations 
I  to  VIII:  A1  (Is),  As  (ls,2s,2p) 

Basis  set  I:  Al(ls,2s,2p,3s), 

As(ls,2s,2p,3s,3p,3d,4s,4p) 

Basis  set  II;  Set  I  plus  Al(3p'^) 

Basis  set  III:  Set  II  plus  As(5s°) 

Basis  set  IV:  Set  III  plus  Al(4s°) 

Basis  set  V:  Set  IV  plus  As(4d°) 

Basis  set  VI:  Set  V  plus  Al(3d°) 

Basis  set  VII:  Set  VI  plus  As(5p°) 

Basis  set  VIII:  Set  VII  plus  Al(4p°) 
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FIG.  1.  Electronic  band  structure  of  AlAs.  The  solid  lines  rep¬ 
resent  the  results  from  Calculation  VII  and  the  dashed  lines  those 
from  Calculation  VIII.  Other  results  in  Tables  I  and  III  and  Figs.  2 
and  3  are  from  Calculation  VII.  The  lattice  constant  is  5.66  A. 

functions  of  the  ionic  states  of  and  As*“  were  obtained 
from  self-consistent,  ab  initio,  atomic  calculations.  The  ra¬ 
dial  parts  of  the  atomic  wave  functions  were  expanded  in 
terms  of  Gaussian  functions.  Sets  of  even-tempered  Gaussian 
exponentials  were  employed  with  a  minimum  of  0.10  for 
both  A1  and  As,  and  a  maximum  of  4.5000  X  10“^  for  A1  and 
1.5000  X  10^  for  As,  in  atomic  units.  We  used  19,  19,  and  17 
orbitals  for  the  s,  p,  and  d  states  for  both  A1  and  As.  A  mesh 
of  60  k  points,  with  proper  weights  in  the  irreducible  Bril- 
louin  zone,  was  employed  in  the  self-consistent  iterations. 
The  computational  error  for  the  valence  charge  was  about 
0.000129  for  34  electrons.  The  self-consistent  potentials  con¬ 
verged  to  a  difference  of  10"^  after  about  30  iterations.  The 
total  number  of  iterations  varied  with  the  input  potentials. 
The  basis  sets  for  each  of  the  self-consistent  calculations 
used  in  the  BZW  procedure  in  this  work  are  listed  in  Table  II. 
The  optimal  basis  set  for  the  AlAs  ealculation  is  the  basis  set 
for  calculation  VII. 


III.  RESULTS 

Figure  1  shows  the  band  structure  of  AlAs  at  high  sym¬ 
metry  points  and  along  high  symmetry  lines  in  the  Brillouin 
zone.  The  calculated  conduction  band  minimum  is  at  the  X 
point  while  the  valence  band  maximum  is  at  F,  in  agreement 
with  experiment.  The  energy  levels  at  high  symmetry  points 


TABLE  III.  Electronic  energies  (eV)  at  high  symmetry  points, 
with  the  top  of  the  valence  band  set  to  0  eV.  Data  are  from  Calcu¬ 
lation  VII  for  a  lattice  constant  of  5.66  A. 
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FIG.  2.  The  calculated  density  of  states  for  AlAs,  as  obtained 
from  the  bands  from  Calculation  VII,  for  a  lattice  constant  of 
5.66  A.  The  inset  shows  the  practically  measurable  band  gap  from 
2.14  to  2.4  eV. 

in  the  Brillouin  zone  are  listed  in  Table  III.  Our  calculated, 
indirect,  fundamental  band  gap  value  of  2.15  eV  agrees  very 
well  with  the  experimental  values  of  2.15  and  2.16  eV,  at 
room  temperature,  and  is  lower,  by  less  than  5%,  than  the 
zero  temperature  values  of  2.23  and  2.25  eV,  as  shown  in 
Table  I.  Similarly,  the  calculated  indirect  band  gap  at  the  L 
point,  2.38  eV,  agrees  very  well  with  the  experimental  find¬ 
ing  of  2.36  eV.  However,  our  direct  band  gap  of  2.35  eV  is 
about  21%  and  25%  lower  than  the  measured  values,  as  dis¬ 
cussed  further  below. 

Figures  2  and  3  show  the  total  (DOS)  and  partial  densities 
of  states  (PDOS)  for  AlAs,  respectively,  as  obtained  from 
Calculation  VII.  The  total  DOS  curves,  particularly  the  inset, 
show  that  the  practically  measurable  band  gap  could  be  any¬ 
where  between  2.14  and  2.4  eV,  where  the  latter  corresponds 
to  a  value  of  0.1  for  the  DOS,  in  excellent  agreement  with 
experiments.  The  lowest  peak  between  -11.8  and 
-9.7  eV  arises  from  the  Al(3p)  and  As(4^)  orbitals.  The 
peaks  appearing  between  -5.6  and  0  eV  stem  from  A1(3j), 
Al(3p),  and  As(4p)  orbitals.  A  significant  hybridization  of 
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FIG.  3.  The  partial  density  of  states  (PDOS)  for  AlAs  from  the 
contribution  of  the  s,  p,  and  d  states  of  A1  and  As  atoms,  respec¬ 
tively.  These  PDOS  are  derived  from  the  bands  from  Calculation 
VII,  for  a  lattice  constant  of  5.66  A. 
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Al(4/?),  Al(3d),  and  As{4d)  orbitals  is  apparent  in  the  PDOS 
curve  in  Fig.  3. 

Our  calculated  electron  effective  mass  of  A1  As,  m*,  at  the 
bottom  of  the  conduction  band  at  the  F  point,  is  0.15  m^, 
where  Wq  is  the  free-electron  mass.  This  result  is  also  in 
excellent  agreement  with  the  experimental  value  of  0.15  Wq 
in  Ref.  32  and  is  21%  higher  than  the  0.124  niQ  estimated 
from  a  fit  to  absorption  data.^^ 

IV.  DISCUSSIONS 

In  comparison  to  previous  density  functional  and  other 
calculations,  this  work  has  resolved  most  of  the  discrepan¬ 
cies  between  experiment  and  theory.  Previous  density  func¬ 
tional  calculations,  as  per  the  content  of  Table  I,  woefully 
underestimated  the  indirect  band  gaps  and  the  direct  one.  The 
GW  calculations  improved  the  agreement  between  theory 
and  experiment  for  the  fundamental,  indirect  gap,  by  reduc¬ 
ing  the  underestimation  to  about  3%  from  the  room  tempera¬ 
ture  value.  For  the  larger  indirect  gap  at  L,  however, "they 
overestimated  significantly  (by  18%  and  28%).  While 
GW-PP'-^  underestimates  the  room  temperature  direct  gap  by 
close  to  8%,  the  GW-QP®  overestimates  it  by  9%.  Vurgaft- 
man  et  al^  underscored  the  difficulties  associated  with  the 
experimental  determination  of  the  direct  gap  at  F.  These  dif¬ 
ficulties  are  apparent  in  the  work  of  Moncmar'  where  high 
impurity  concentrations  caused  a  broadening  of  exciton  spec¬ 
tra.  Consequently,  the  well-defined  dip  in  the  spectrum,  that 
would  have  indicated  the  direct  gap,  is  replaced  by  a  much 
broader  dip  from  the  uncertain  location  of  which  one  extracts 
the  value  of  the  direct  gap.  We  expect  the  uncertainties  asso¬ 
ciated  with  the  determination  of  this  direct  gap  to  contribute 
to  the  unusual  deviation  of  our  calculated  value  from  experi¬ 
ment.  Further,  given  that  this  gap  is  not  the  minimum  one, 
GGA-BZW  calculations  are  not  necessarily  expected  to  de¬ 
scribe  it  correctly  from  a  ground  state  theory  as  it  does  for 
the  fundamental  bang  gap.^°“^^ 

A  central  issue  in  these  discussions  relates  to  the  physical 
and  mathematical  reasons  that  explain  the  success  of  the 
BZW  method  where  many  others  have  failed,  including 
some  rather  sophisticated  approaches  that  go  beyond  density 
functional  theory.  The  answer  resides  in  a  thorough  under¬ 
standing  of  basis  set  related  effects  on  the  outcomes  of  these 
calculations.  Specifically,  as  discussed  by  Bagayoko  et  al,^^ 
LCAO  type  calculations  have  to  ascertain  the  completeness 
of  the  basis  set.  In  doing  so,  one  typically  includes  diffuse 
orbitals,  with  relatively  small  exponents,  to  accommodate  the 
charge  redistribution  in  the  molecular  or  solid  state  environ¬ 
ments  as  compared  to  those  for  atoms  or  ions.  Additionally, 
one  has  to  provide  for  angular  symmetries  that  may  arise  in 
molecular  or  solid  environments.  Most  works  cited  herein 
appear  to  have  handled  these  two  points  correctly. 

There  is  still  a  nontrivial  basis  set  related  effect  distinct 
from  the  issues  of  angular  symmetry  and  of  diffuse  orbitals. 
The  preoccupation  for  guaranteeing  the  completeness  of  the 
basis  set  explains  the  use  of  as  large  basis  sets  as  possible,  as 
long  as  linear  dependency  is  avoided.  Consequently,  no  par¬ 
ticular  limit  is  placed  on  the  size  of  the  basis  set.  This  size 
often  varies  widely  from  one  group  of  investigators  to  an¬ 


other,  even  for  the  same  system  and  for  similar  orbitals  (i.e., 
Gaussians).  The  Rayleigh  theorem,^®’^^  however,  indicates 
that  serious  problems  could  be  associated  with  arbitrarily 
large  basis  sets.  The  statement  of  this  rigorous,  mathematical 
theorem  follows:  Let  an  eigenvalue  equation  be  solved  twice, 
by  an  LCAO  method,  respectively,  with  N  and  (N+l) 
orbitals — such  that  the  (N+l)  orbitals  of  Calculation  II  in¬ 
clude  all  the  N  orbitals  of  Calculation  I  plus  an  additional 
one — and  let  the  eigenvalues  from  the  two  calculations  be 
ordered  from  the  lowest  (sf  and  I  and  II,  respec¬ 
tively)  to  the  highest  (ej^  and  respectively),  then  the 

theorem  states  that  for  all  i^N.  Alternatively, 

this  theorem  states  that  a  variational  eigenvalue,  upon  an 
increase  of  the  basis  set  (and  hence  of  the  dimension  of  the 
matrix),  is  not  increased.  It  either  stays  the  same  (if  it  is 
equal  to  the  exact,  physical  eigenvalue  of  the  matrix)  or  it 
decreases  to  approach  this  exact  eigenvalue  from  above.  This 
theorem  therefore  points  to  the  need  to  utilize  as  large  a  basis 
set  as  possible  in  order  to  lower  variational  eigenvalues  to 
reach  the  exact  ones.  Such  a  large  basis  set  is  also  needed  to 
ensure  completeness.  Bagayoko,  Zhao,  and  Williams^®'^* 
identified  a  basis  set  and  variational  effect  stemming  from 
any  lowering  of  unoccupied  energy  levels  or  bands,  for  mol¬ 
ecules  or  solids,  respeetively,  when  the  occupied  levels  or 
bands  have  converged  (vis-a-vis  an  increase  of  the  basis  set). 
The  effect  can  exist  in  any  variational  calculations  for  which 
the  occupied  and  unoceupied  states  are  not  treated  in  identi¬ 
cal  manners.  Specifically,  for  all  variational  calculations  of 
which  we  know,  the  occupied  and  unoccupied  states  are 
treated  in  fundamentally  different  manners:  in  the  iterative 
process  toward  self-consistency,  the  eharge  density  is  con¬ 
structed  or  reconstructed  using  only  the  wave  functions  of  the 
occupied  states.  Hence,  only  the  wave  functions  of  these 
states  affect  the  charge  density,  the  potential,  and  the  Hamil¬ 
tonian  itself,  i.e.,  the  physics  of  the  system.  It  is  therefore 
necessary  to  guarantee  that  an  applicable  basis  set  is  large 
enough  to  ensure  the  convergence  of  the  occupied  levels  or 
bands  vis-a-vis  the  size  of  the  basis  set.  To  do  so,  Bagayoko, 
Zhao,  and  Williams  (BZW)  introduced  their  method  that  re¬ 
quires  the  utilization  of  a  minimum  basis  set  for  the  first, 
self-consistent  calculation.  The  outcomes  of  this  first  calcu¬ 
lation  are  generally  not  the  correct  solutions  for  the  physical 
system  under  study.  The  method  therefore  requires  subse¬ 
quent  calculations  where  the  basis  set  is  methodically  aug¬ 
mented  as  described  above  and  elsewhere.^'*’^^  The  occupied 
energy  levels  or  bands  of  a  calculation  are  compared  to  those 
of  the  previous  one  that  has  a  smaller  basis  set.  This  proeess 
continues  until  calculation  N  and  (A+1)  have  the  same  oc¬ 
cupied  energy  levels  or  bands,  within  applicable  uncertain¬ 
ties.  We  then  select  the  outputs  of  calculation  N  as  the  physi¬ 
cal  description  of  the  system.  Indeed,  in  light  of  the  Rayleigh 
theorem,  some  unoccupied  energy  levels  or  bands  of  Calcu¬ 
lation  {N+ 1)  are  generally  lower  as  compared  to  their  coun¬ 
terpart  from  Calculation  N.  The  BZW  method  ascribes  any 
such  additional  lowering  of  an  unoccupied  band,  while  the 
occupied  ones  do  not  change,  to  mathematical  artifacts  stem¬ 
ming  directly  from  the  Rayleigh  theorem.  While  the  method 
applies  to  all  variational  calculations  of  the  Rayleigh-Ritz 
type  that  employ  an  LCAO  approach,  it  is  particularly  perti- 
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nent  for  density  functional  calculations.  Indeed,  density  func¬ 
tional  theory  is  fundamentally  a  ground  state  theory.  The 
BZW  method  verifiably  guarantees  the  proper  description  of 
the  ground  state  with  a  basis  that  is  sufficiently  large.  Varia¬ 
tional  calculations  that  do  not  employ  the  method  generally 
have  large  enough  basis  sets  to  describe  the  applicable  sys¬ 
tem.  However,  they  do  not  actually  verify  that  the  size  of  the 
basis  set  converged  vis-a-vis  the  description  of  the  occupied 
states  and  they  do  not  avoid  Rayleigh-theorem  related  low¬ 
ering  of  some  unoccupied  levels  or  bands  on  purely  math¬ 
ematical  ground — given  that  the  Hamiltonian  (i.e.,  the  phys¬ 
ics)  does  not  change  once  the  optimal  basis  size  of  the  BZW 
method  is  attained. 

V.  CONCLUSION 

In  conclusion,  we  successfully  performed  ab  initio,  self- 
consistent  GGA-BZW  calculations  for  the  band  structure  of  a 


zinc-blende  AlAs  semiconductor.  Our  calculated  band  gaps, 
except  for  the  direct  gap  at  F,  are  in  excellent  agreement 
with  measured  values.  Our  calculated  effective  mass  also 
agrees  very  well  with  a  directly  determined  experimental 
value  and  disagrees  with  an  estimate  derived  from  fitting  to 
absorption  data. 
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We  report  structural,  elastic,  and  electronic  properties  of  selected,  deformed,  single-wall  carbon  nanotubes 
under  uniaxial  strain.  We  utilized  a  generalized  gradient  approximation  potential  of  density  functional  theory 
and  the  linear  combination  of  atomic  orbital  formalism.  We  discuss  bond-lengths,  tubule  radii,  and  the  band 
gaps  as  functions  of  tension  and  compression  strain  for  carbon  nanotubes  (10,  0),  (8,  4),  and  (10,  10)  which 
have  chiral  angles  of  0,  19.1,  and  30  deg  relative  to  the  zigzag  direction.  We  also  calculated  the  Young’s 
modulus  and  the  in-plane  stiffness  for  each  of  these  three  nanotubes  as  representatives  of  zigzag,  chiral,  and 
armchair  nanotubes,  respectively.  We  found  that  these  carbon  nanotubes  have  unique  structural  properties 
consisting  of  a  strong  tendency  to  retain  their  tubule  radii  under  large  tension  and  compression  strains. 

DOI:  10.1103/PhysRevB.71.205410  PACS  number(s):  73.22. -f,  6 1.46.4- w 


Single-wall  carbon  nanotubes  (SWCNTs)  can  be  viewed 
as  rolled-up  graphene  sheets  that  have  a  diameter  at  the  order 
of  1  nm.  They  have  properties  such  as  high  current  density, 
high  elasticity,  and  stifftiess  unparalleled  by  other  materials. 
Their  potential  applications  range  from  that  in  building  sky¬ 
scrapers  and  elevator  cables  to  the  ones  in  very  tiny  electrical 
circuits  and  machines.’  These  materials,  however,  are  too 
small  for  many  conventional  measurements.  This  situation 
underscores  the  possible  importance  of  theoretical  studies, 
including  the  one  reported  here  that  focuses  on  structural, 
elastic,  and  electronic  properties  of  selected  single-wall  car¬ 
bon  nanotubes  under  uniaxial  strain. 

In  the  last  several  years,  tight-binding  calculations  have 
been  extensively  used  to  study  the  structural  and  electronic 
suaicture  of  carbon  nanotubes.^"^  Tight-binding  approxima¬ 
tions  based  on  the  symmetry  of  the  honeycomb  lattice  of 
graphite  predicted  that  SWCNTs  could  be  semiconducting  or 
metallic  depending  on  their  chirality  (n,m).  The  tight- 
binding  model  has  been  able  to  provide  good  estimates  of  the 
basic  electronic  structure  of  SWCNTs.  However,  curvature- 
related  effects  and  the  hybridization  of  different  electronic 
states  of  graphite  could  lead  to  structural  and  electronic  prop¬ 
erties  that  are  substantially  different  from  the  result  of  tight- 
binding  calculations.  Zigzag  (n,0)  (where  n  is  a  multiple  of 
3)  SWCNTs  which  were  predicted  to  be  metallic  from  tight- 
binding  calculations  were  found  to  have  small  energy  gaps.’ 
Several  theoretical  groups  have  studied  the  elastic  properties 
of  carbon  nanotubes.  Their  approaches  include  simulations 
with  realistic  many-body  potentials,®  the  empirical  force- 
constant  method,®  tight-binding  formalisms,®’’®  pseudopoten¬ 
tial  calculations  with  local  density  approximation 
potentials,”  and  Bom’s  perturbation  technique  within  a  lat¬ 
tice  dynamics  model.” 

The  aim  of  this  work  is  therefore  to  study  the  aforemen¬ 
tioned  properties  of  SWCNTs,  utilizing  ab  initio  quantum 
computations.  We  specifically  report  on  structural,  elastic, 
and  electronic  properties  of  SWCNTs  under  uniaxial  strain. 
Recent  tight-binding  calculations  have  led  to  values  of  bond 
lengths  and  radii,  band  gaps,  and  Fermi  levels  as  functions  of 
strain.®  The  dependence  of  tight-binding  results  upon  the  se¬ 


lected  parameters  partly  suggests  the  present  work  that  em¬ 
ploys  an  ab  initio  approach. 

■Three  main  features  characterize  our  first-principles  com¬ 
putational  method.  The  first  of  these  features  is  the  choice  of 
tlie  potential.  We  utilized  the  generalized  gradient  approxi¬ 
mation  (GGA)  potential  of  Perdew  and  Wang.’®  This  density 
functional’'*"’®  potential  goes  beyond  the  local  density  ap¬ 
proximation  (LDA).  We  also  performed  LDA  calculations 
for  carbon  nanotube  (10,  0).  The  second  feature  of  our 
method  stems  from  employing  the  linear  combination  of 
atomic  orbital  (LCAO).  The  third  and  distinctive  feature  of 
our  work  resides  in  our  use  of  optimal  basis  sets  as  per  the 
Bagayoko,  Zhao,  and  Williams  (BZW)  procedure.”"®  As 
explained  elsewhere,  this  procedure  avoids  a  basis  set  and 
variational  effect  inherently  associated  with  variational  cal¬ 
culations  that  employ  a  basis  set  and  leads  to  the  calculated 
band  gaps  in  very  good  agreement  with  experiments.”"® 
With  tlie  above  method,  we  solved  the  Kohn-Sham’^’’®  equa¬ 
tion  self-consistently.  Self-consistency  was  followed  by  the 
calculations  of  the  total  energies.  Details  of  these  steps,  in¬ 
cluding  the  Kohn-Sham  equation  and  the  expression  for  the 
total  energy,  are  fully  described  in  the  literature.’^"® 

In  tlie  LCAO  method,  we  expanded  the  electronic  eigen¬ 
function  of  the  many-atom  system  as  a  linear  combina¬ 
tion  of  atomic  wave  functions.”  These  input  functions  result 
from  ab  initio  calculations  for  atomic  species  that  are  present 
in  the  system.  For  the  calculations,  the  C(1j)  state  was  used 
as  the  core  state.  For  the  (10,  0)  and  the  (8,  4)  tubes,  the 
C(2s2p)  states  were  used  as  filled  and  partially  filled  valence 
states.  The  C(3s3p)  orbitals  were  the  unfilled  electron  states 
that  were  used  to  augment  the  basis  set  for  the  calculations. 
For  the  (10,  10)  tube,  the  C{2s2p)  states  were  used  as  filled 
and  partially  filled  valence  states  and  the  C(35)  as  the  un¬ 
filled  state.  The  empty  C(3p)  state  was  dropped,  in  the  case 
of  (10,  10),  due  to  convergence  difficulties. 

Uniaxial  strain  was  simulated  by  linearly  scaling  the  atom 
positions  along  the  tube  axis  in  the  carbon  nanotube.  To  find 
how  the  tube  radius  changes  with  uniaxial  strain,  after  the 
tube  axis  was  scaled,  the  radius  was  identified  from  the  total 
energy  minimization  procedure.  Namely,  for  a  given  nano- 
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TABLE  II.  Coefficients  m  of  linear  variation  of  normalized  bond 
lengths  with  compression  (negative)  and  tensile  (positive)  strain 
(i.e.,  e  expressed  as  -0.04  and  +0.06  for  4%  and  6%  compression 
and  tensile  strain,  respectively).  Al/lo=me,  where  1  is  the  a,  b,  or 
c  bond  length.  The  Young’s  modulus  Y  is  in  TPa.  The  in-plane 
stiffness  C  is  in  J/m^. 


Nanotubes 

m  for  bond 

a 

m  for  bond 
b 

m  for  bond 

c 

Y 

(TPa) 

C 

(J/m®) 

(10,0) 

0.252 

0.252 

1,000 

1.47 

340 

(8.4) 

0.036 

0.572 

0.893 

1.10 

267 

(10,10) 

0.746 

0.746 

0“ 

0.726 

272 

“Bond  c,  for  (10,10),  is  along  the  circumference  of  the  tube. 


Oq  and  bo  and  1.406  A  for  Cq,  that  is  about  1%  smaller  than 
the  results  of  tlie  ab  initio  GGA  calculations.  The  ab  initio 
LDA  calculations  for  SWCNT  (10,  0)  also  found  tliat  the 
radii  of  the  nanotube  did  not  change  under  a  substantially 
large  uniaxial  strain.  Among  the  previous  ab  initio  calcula¬ 
tions  of  the  elastic  properties  of  carbon  nanotubes,  Sanchez- 
Portal  et  al.  utilized  a  minimal  basis  set  of  one  s  and  three  p 
orbitals  per  carbon  atom  and  performed  LDA  calculations.^ 
In  tlieir  calculations,  they  used  pseudoatomic  orbitals.^  Van 
Lier  et  al.  utilized  ab  initio  Hartree-Fock  6-3 IG  method  and 
closed  nanotube  models  in  tlteir  simulations Sanchez- 
Portal  et  al.^  and  Van  Lier  et  al.~^  reported  relatively  small 
values  (from  0.14  to  0.19)  of  the  Poisson  ratio  for  their  cal¬ 
culated  carbon  nanotubes,  which  also  indicated  the  diameter 
rigidity  of  the  carbon  nanotubes.  These  previous  ab  initio 
calculations  utilized  different  computational  methods,  such 
as  a  minimal  basis  set  and  tlie  Hartree-Fock  method,  and 
their  results  are  slightly  different  from  ours.  As  demonstrated 
in  one  of  our  previous  publications,^®  the  use  of  a  minimal 
basis  set  may  not  be  sufficient  to  obtain  a  highly  accurate 
solution  of  the  calculated  electronic  stmcture  of  the  carbon 
nanotubes,  which  could  explain  in  part  the  difference  be¬ 
tween  our  results  and  those  of  Sanchez-Portal  et  al  The  fun¬ 
damental  differences  between  the  Hartree-Fock  method  and 
density  functional  methods,  i.e.,  the  inclusion  of  the  elec¬ 
tronic  correlation  effects  in  the  latter,  partly  explain  the  dif¬ 
ference  between  our  results  and  those  of  Van  Lier  et  al  In 
our  ab  initio  calculations,  we  utilized  extended  atomic  basis 
sets  and  performed  both  GGA  and  LDA  computations. 

The  calculated  Young’s  moduli  Y  (in  TPa)  for  the  tubes 
are  shown  in  Table  H.  We  recall  that  the  radius  included  in 
the  formula  for  the  Young’s  modulus  is  from  the  center  of  the 
tube  to  the  outer  circumference.  This  value  of  tlie  radius  led 
to  Young’s  moduli  close  to  the  experimentally  found  value  of 
approximately  1  TPa.®  The  (10,  0)  and  (10,  10)  tubes  have 
the  highest  and  lowest  Young’s  moduli,  respectively.  This 
trend  is  expected  due  to  the  bond  geometry  of  the  tubes 
according  to  chirality.  All  three  bonds  in  the  (10,  0)  tube 
have  a  significant  component  parallel  to  the  axis  along  which 
strain  is  applied.  In  the  case  of  (10,  0),  the  c  bond  is  entirely 
along  the  axis.  In  contrast,  one  of  the  bonds  of  the  (10,  10) 
tube  is  entirely  along  the  circumference;  it  is  not  expected  to 
oppose  any  resistance  neither  to  compression  nor  to  tensile 
strain.  For  this  reason,  nanotube  (10,  10)  is  intuitively  the 


e  (%) 

FIG.  3.  Band  gaps  of  indicated  carbon  nanotubes  as  functions  of 
compression  and  tensile  strain. 

easiest  to  stretch  or  to  compress,  as  confirmed  by  our  results 
in  Table  H. 

We  also  calculated  the  in-plane  stiffness  C  for  the  carbon 
nanotubes.  It  is  340,  267,  and  272  J/m^  for  (10,  0),  (8,  4), 
and  (10,  10),  respectively.  The  in-plane  stiffness  of  these 
nanotubes  exhibited  a  dependence  on  their  chirality,  when 
the  nanotubes  are  strained  along  the  mbe  axis  (z-direction). 
Xiao  and  Liao  reported  an  average  in-plane  stiffness  C  of 
328  J/m^  for  graphene,  using  the  second-generation  Brenner 
potential  in  their  simulation  of  carbon  nanotubes 

The  calculated  electronic  properties  for  the  unstrained 
nanotubes  basically  reproduced  the  results  of  Zhao  et  aO^ 
for  (10,  0)  and  (8,  4),  using  the  optimal  basis  sets  of  the 
BZW  method,  which  gave  converged  results  for  the  calcu¬ 
lated  electronic  structure  that  included  the  occupied  electron 
states  as  well  as  the  unoccupied  ones  near  the  Fermi  level. 
These  tubes  were  found  to  be  semiconductors.  These  authors 
provided  plots  of  the  band  structures  of  these  tubules.^®  The 
electron  energy  bands  of  carbon  nanotube  (10,  10)  present  a 
semimetallic  property  as  reported  in  previous  publications. 
Figure  3  shows  the  nonlinear  variation  of  band  gaps  of  (10, 
0)  and  (8,  4)  with  strain.  For  both  nanotubes,  the  band  gap 
decreases  with  compression  strain  and  increases  with  tensile 
strain,  for  strain  values  smaller  than  or  equal  to  4%.  While 
the  band  gap  for  (8, 4)  reaches  a  minimum  for  a  compression 
strain  of  4%,  that  for  (10,  0)  exhibits  a  maximum  for  a  ten¬ 
sile  strain  of  4%.  Our  results  for  the  variation  of  the  band  gap 
with  strain  qualitatively  agree  with  the  findings  from  tight 
binding^’  but  are  quantitatively  different,®  particularly  in  the 
magnitudes  of  the  band  gaps.  As  apparent  in  Fig.  3,  our 
results  indicate  that  nanombe  (10,  0)  become  metallic  at  6% 
compression  strain.  It  was  noticed  that  the  carbon  nanotubes 
may  collapse  for  large  strains,®^’®®  a  process  that  we  cannot 
simulate  at  present  using  the  ab  initio  quantum  calculations, 
because  the  required  computations  are  beyond  our  current 
computation  capability.  Therefore,  the  obtained  band  gaps 
for  large  strains  may  be  a  theoretical  simulation.  A  priori,  the 
general  behavior  of  these  band  gaps  with  strain,  particularly 
from  values  below  4%,  is  expected  on  the  basis  of  the  in¬ 
crease  or  decrease  of  the  overlap  between  atomic  sites,  for 
compression  or  tensile  strains,  respectively.  The  consequent 
broadening  or  flattening  of  the  bands,  respectively,  affects 
the  gaps. 

A  discussion  of  our  results  for  the  elastic  properties  is 
partly  limited  by  the  dearth  of  experimental  data  and  by  the 
rather  large  uncertainty  associated  with  currently  available 
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Abstract:  We  performed  ab-initio  density  functional  simulations  to  study  the  structural  and 
growth  properties  of  short  carbon  nanobells.  We  used  a  real  space  approach  and  the  linear 
combination  of  atomic  orbitals  (LCAO)  formalism.  In  the  nitrogen-doped  carbon  nanobells,  the 
nitrogen  atoms  that  are  attracted  to  the  open-edge  sites  of  the  carbon  nanobells  play  an  important 
role  in  the  growth  of  the  short  carbon  nanostructures.  We  present  the  calculated  electronic 
stracture  of  the  short  nanobells.  The  calculated  local  density  of  states  of  the  nanobells  revealed 
field  emission  characteristics  that  agi'ee  with  experimental  observations. 


I.  Introduction 

The  cai'bon  nanostructures  present  interesting 
properties  with  a  great  potential  of 
applications.  Particularly,  short  carbon 
nanotubes  or  carbon  nanobells  present 
excellent  field-emission  properties  that  have 
attracted  a  great  interest  for  applications.^'^ 
Experiments  demonstrated  that  the  nitrogen- 
assisted  synthesis  can  grow  carbon  nanofibers 
on  a  large  scale.  '  Such  carbon  nanofibers 
exhibit  a  "bamboo-like"  structure. 
Distinctively,  a  great  part  of  the  bamboo-like 
nanofibers  consists  of  short  carbon  nanobells. 
Individual  nanobeils  are  self-contained  and 
stack  one  on  top  of  the  other  to  create  a  long 
nanofiber.  The  carbon  nanobells  may  be 
viewed  as  short  nanotubes  such  that  their 
lengths  are  of  the  same  order  as  their 
diameters.  The  electronic  structure  of  the 
short  carbon  nanobells  is  substantially 
different  from  that  of  the  pure  and  long 
carbon  nanotubes  because  of  the  quantum 
effects  inherent  to  their  sizes.  The  nanobeils 
exhibit  novel  electron  field  emission 


properties  with  a  tum-on  field  of  electron 

Q  1  <5 

emission  as  low  as  0.8E//^?r.’  We 
performed  ab-initio  density  functional 
calculations,  aiming  to  understand  the 
growth,  structural,  and  electronic  properties 
of  short  carbon  nanobeils. 

II.  Method 

Although  the  experimental  synthesis  of 
carbon  nanobeils  produced  samples  of 
various  sizes,  measurements  revealed  that  the 
atomic  structure  of  the  bamboo-like 
morphology  in  thick  nanobeils  is  similar  to 
that  of  the  thinner  ones.®  Experiments  also 
revealed  that  adjacent  nanobeils  do  not  have  a 
firm  contact,  but  instead  the  closed  end  of  one 
nanobell  is  weakly  inserted  into  the  open  end 
of  another.  Single  nanobeils  can  be  easily 
separated  from  other  part  of  the  nanofiber. 
Hence,  our  ab-initio  calculations  will  focus 
on  single  nanobeils.  The  weak  interaction 
between  adjacent  carbon  nanobeils  may  be 
included  as  a  perturbation  in  further  studies  in 
modeling  the  bamboo-like  morphology  of  the  • 
nanofibers.  We  constructed  a  prototypical 
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in  the  inner  rings,  away  from  the  open-edge 
of  the  nanobell.  This  result  indicates  that  the 
effect  of  further  relaxation  of  the  third  atomic 
ring  and  inner  rings  would  be  insignificant. 

(B).  A  Possible  Growth  Mechanism  of 
Carbon  Nanobells.  The  growth  mechanism 
of  the  carbon  nanobells  remains  a  difficult 
problem.  There  are  probably  two  major 
growth  models.’  The  first  one  assumes  that 
the  carbon  atoms  are  added  at  the  open-ends 
of  the  nanobells.  The  second  one  involves 
the  C2  absorption  process  that  is  assisted  by 
the  pentagonal  defects  on  the  closed  caps.’ 
Although  both  models  are  very  interesting, 
we  have  technical  difficulties  using  density 
functional  computations  to  simulate  the 
second  growth  model.  We  successfully 
performed  local  density  functional 
computations  to  test  the  first  growth 
mechanism.  Because  of  the  active  dangling 
bonds  of  carbon  atoms  at  the  open-edge  of 
nanobells,  there  is  a  high  possibility  that 
carbon  atoms  can  be  attracted  to  these  sites 
for  growth  to  occur.  We  calculated  the  total 
energies  in  the  following  two  cases.  In  the 
first  case,  a  ring  of  20  carbon  atoms  attaches 
to  the  open-edge  of  the  carbon  nanobell  at  a 

bond-length  of  1.415  A.  In  the  second  case, 
these  20  carbon  atoms  are  free.  The 
calculated  total  energy  in  case  I  is 
substantially  lower  than  the  corresponding 
value  in  case  11.  The  total  energy  difference, 
which  may  also  be  defined  as  the  cohesive 
energy  of  the  carbon  atoms  in  the  first  growth 
model,  is  4.7  eV/atom.  In  this  growth  model, 
see  Fig.  1,  C2  dimers  are  absorbed  at  the 
active  dangling  bond  sites  at  the  open  end  of 
the  carbon  nanobell.  A  C2  dimer  that  deposits 
on  the  open-edge  of  the  bell  forms  one 
covalent  bond  on  the  same  atomic  ring.  Each 
of  the  C  atoms  also  forms  one  covalent  bond 
with  the  C  atom  of  the  next  atomic  ring  and 
remains  an  active  dangling  bond  towards  the 
open-space.  This  grow  model  is  consistent 
with  previous  analysis  of  experimental 


results.”’’^’  Although  we  cannot  rule  out 
the  second  growth  model  (growth  on  the 
closed  cap  of  carbon  nanobells),  our 
calculations  presented  a  clear  evidence  that 
the  first  growth  mechanism  is  highly  possible. 

(C).  Nitrogen-Doped  Carbon  Nanobell. 

Experimental  results  indicated  that  the  growth 
of  the  short  carbon  nanobells  is  highly 
dependent  on  the  nitrogen  concentration  in 
the  gas  mixture  during  the  synthesis.®’’” 
Without  nitrogen  atoms  in  the  growth  gas 
mixture  under  the  same  conditions,  long  and 
pure  carbon  nanotubes,  without  the  “bamboo¬ 
like”  morphology,  are  produced.  However,  it 
was  not  clear  why  the  nitrogen  atoms  could 
turn  the  growth  of  the  would-be  long  carbon 
nanotubes  into  that  of  short  carbon  nanobells. 
We  performed  total  energy  calculations  to 
study  the  nitrogen-doped  carbon  nanobells 
(CNx  nanobells).  Since  nitrogen  atoms  can 
also  form  the  planar  sp^-hybrid,^^  we  studied 
the  substitutional  doping  of  nitrogen  atoms  in 
the  CNx  nanobells.  The  substitutional  doping 
of  nitrogen  atoms  in  CNx  nanobells  was  also 
proposed  from  the  analysis  of  experimental 
results.’^ 

We  first  replaced  carbon  atoms  with 
nitrogen  atoms  at  the  open-edge  of  the 
prototypical  nanobell.  We  compared  the  total 
energies  in  the  two  cases.  In  the  first  case, 
computations  are  carried  out  for  the  pristine 
carbon  nanobell  in  the  presence  of  ten  free 
atoms  of  nitrogen.  In  the  second  case,  ten 
nitrogen  atoms  replace  ten  carbon  atoms  on 
the  open-edge  and  the  substitutionally  N- 
doped  nanobell  is  in  the  presence  of  ten  free 
atoms  of  carbon.  The  calculated  total  energy 
in  the  second  case  is  lower  than  that  in  the 
first  case  by  an  amount  of  0.5  eV  per  atom, 
without  relaxing  the  structure.  We  then 
performed  total  energy  calculations  and 
identified  the  atomic  positions  of  the  relaxed 
stracture  of  the  nitrogen-doped  nanobell  from 
the  total  energy  .minimization.  The  calculated 
C-N  bond  length  in  the  first  atomic  ring  at  the 
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about  0.08  states  per  eV  per  atom.  Thus, 
electrons  will  be  emitted  from  the  atomic 
sites  near  the  open-edge  of  the  nanobell  in 
field-emission  experiments.  The  LDOS  of 
Fig.  2  exhibit  a  nearly  metallic  behavior.  The 
real  nanobell  samples  in  experimental  studies 
mostly  involve  multilayers  of  graphite  sheets 
and  the  nanobells  stacked  one  on  top  of  the 
other.  The  weak  interactions  between  graphite 
layers  in  these  nanobells  and  that  between  the 
bells  may  further  broaden  the  electronic 
energy  levels.  Consequently,  a  metallic 
behavior  may  be  observable  in  the 
measurements  of  the  electronic  structure  of 
N-doped  nanobells. 

IV.  Conclusion 

In  conclusion,  we  performed  ab-initio 
density  functional  calculations  to  simulate  the 
growth,  structural,  and  electronic  properties 
of  prototypical  carbon  nanobells.  In  the 
nitrogen-doped  carbon  nanobells,  nitrogen 
atoms  that  are  attracted  to  the  open-edge  of 
the  nanobells  may  play  a  role  of  stopper  of 
the  growth  of  the  nanostructures.  The 
calculated  local  densities  of  states  of  the  CNx 
nanobells  indicate  that  electrons  are  most 
likely  emitted  from  the  atomic  sites  near  the 
open-edge  of  the  carbon  nanobell  in  field- 
emission  experiments.  This  result  agrees  with 
experimental  observations. 
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Abstract.  We  present  a  universal  relation  between  the  densities  of  states  near  Van  Hove 
singularities  and  the  effective  electron  masses  in  1-dimensional  (1-D)  semiconductors. 
The  relation  can  be  utilized  as  a  new  method  to  determine  the  effective  masses  of  charge 
carriers  of  1-D  semiconductors  in  theoretical  calculations  and  in  experimental 
measurements.  The  calculated  results,  utilizing  the  relation  for  1-D  single-walled  carbon 
nanotubes  (SWCNTs),  agree  well  with  that  of  the  conventional  calculations  when  both 
approaches  utilize  the  outputs  of  ab-initio  density  functional  computations. 


1.  Introduction. 

One-dimensional  (1-D)  materials,  such 
as  carbon  nanotubes  (CNTs),  have 
emerged  as  attractive  materials  for 
applications  in  molecular-scale 
electronics.*'"  The  electronic  properties 
of  1-D  CNTs  have  been  extensively 
studied  in  recent  years  using  various 
experimental  tools  and  theoretical 
calculations."’*"’"'  Particularly,  scanning 
tunneling  microscopy  (STM)  is  a 
favorable  tool  for  measuring  atomic 
structures,  density  of  states  (DOS)  of 
electrons,  Van  Hove  singularities,  and 
energy  band  gaps.  However,  there  are 
some  other  physical  quantities  that 
require  further  studies.  One  of  them  is 
the  effective  mass  of  charge  carriers 
(electrons  or  holes),  an  important 
physical  quantity  that  characterizes  the 
electronic  transport  properties  of 
materials.  Here,  we  present  a  universal 
relation  between  the  densities  of  states 
of  electrons  near  Van  Hove  singularities 
and  effective  electron  masses  in  1-D 


materials.  This  relation  can  be  used  as  a 
new  method  to  determine  the  effective 
masses  of  charge  carriers  in  any  1-D 
semiconductors. 


II.  Densities  of  States  and  Effective 
Masses  in  1-D  Materials. 

The  electron  density  of  states  (DOS)  per 
unit  cell,  G(e),  of  a  1-D  material  can  be 
calculated  by  the  following  expression.'' 


* 


ds.. 


dk 


(1) 


where  is  the  energy  of  the  i-th 

electron  band  at  point  k  in  the 
irreducible  Brillouin  zone;  and  L  is  the 
length  of  the  unit-cell. 

At  an  extremal  point  of  an 
electron  energy  band,  the  first  derivative 
of  the  band  energy  with  respect  to  the 
wave  vector  k  is  zero,  that  is  /dk  = 

0.  A  Van  Hove  singularity  occurs  at 
this  point..  We  can  then  approximate  the 
band  dispersion  as  a  parabolic 
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which  the  band  dispersion  follows  Eq. 
(5).  The  density  of  states  and  the 
effective  mass  of  this  single  band  has 
been  discussed  in  the  previous  section. 
We  then  identify  the  relation  as, 

^2(g  +  gi)  .g. 

2D,(£  +  s,)j 


where  and  D2(£)  are  the  average 
effective  mass  and  the  density  of  states 
of  the  doubly  degenerate  bands, 
respectively. 


(10) 


where  /ni''  = - and  = - . 

The  above  discussions  can  also 
be  extended  to  three  or  more  energy 
bands  that  are  degenerate  at  an  extremal 
point,  but  are  not  degenerate  elsewhere. 


IV.  Applications. 

Although  there  is  a  wealth  of 
experimental  results  regarding  the 
density  of  states  of  electrons  in  single- 
walled  carbon  nanotubes  (SWCNTs),”’“‘ 
the  study  of  the  effective  masses  of  the 
charge  carriers  in  these  materials  is  far 
from  complete.  As  a  test  case  of  the  new 
method  to  study  the  effective  masses  of 
the  charge  carriers  in  1-D  materials^  we 
applied  Eq.(7)  or  (9)  to  some  SWCNTs. 
We  list  some  of  the  results  in  Table  1, 
where  mo  is  the  free  electron  mass. 
in  Table  1  refers  to  the  first  valence  band 
near  the  Fermi  level  (Ep),  whereas 
refers  to  the  first  conduction  bands  from 
Ep.  Em  is  the  energy  (in  eV)  of  the  top 
of  a  valence  band  or  the  bottom  of  a 
conduction  band  near  its  Van  Hove 


singularity.  The  energy  of  the  top  of  the 
first  valence  band  has  been 

chosen  to  be  zero  in  Table  1.  "Deg" 
refers  to  the  degeneracy  of  the  bands. 
Dm  (in  states/eV)  is  the  peak  value  of 
the  density  of  states  at  the  Van  Hove 
singularity;  it  was  used  to  calculate  the 
effective  mass  m*  using  the  universal 
relation,  Eq.(7)  or  (9).  We  utilized  the 
effective  mass  of  charge  carriers  of 
SWCNT  (10,  0)  at  as  the  single 
reference  for  the  calculations  of  the 
effective  masses  (m*)  of  other 
SWCNTs.  As  a  comparison,  we  also 
calculated  the  effective  mass  m*  using 
the  results  of  the  conventional  method 
for  which  we  utilized  the  ab-initio 
electronic  structure  calculations.  We 
performed  self-consistent,  ab-  initio 
density  functional  calculations  to  obtain 
the  electronic  structure  of  the  SWCNTs, 
using  a  linear  combination  of  atomic 
orbital  (LCAO)  method.^'’®  Particularly, 
the  ab-initio  calculations  for  the 
electronic  structure  of  nonsymmorphic 
SWCNTs  (8,  4)  and  (10,  5)  are  very 
difficult  tasks,  since  there  are  112  and 
140  atoms  in  their  tubule  unit  cells, 
respectively.  The  calculated  results  of 
the  effective  masses  7?7*of  the  charge 
carriers  in  SWCNTs  utilizing  the 
universal  relation  Eq.  (7)  or  (9)  agree 
very  well  with  those  from  the 
conventional  method  of  ab-initio 
electronic  structure  calculations,  as 
demonstrated  in  Table  1. 
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Abstract 

We  report  calculated,  electronic  properties  of  wurtzite  zinc  oxide  (w-ZnO).  Unlike  many 
previous  theoretical  works,  our  linear  combination  of  atomic  orbital  (LCAO) 
calculations,  implemented  following  the  Bagayoko,  Zhao,  and  Williams  (BZW)  method, 
employed  a  local  density  approximation  (LDA)  potential  to  obtain  band  gaps  and  other 
results  in  agreement  with  experiments.  We  discuss  the  band  structures,  the  calculated 
direct  band  gap  (3.2  to  3.3  eV),  and  the  electron  effective  mass  (0.25  to  0.26  mo). 

PACS  Numbers  71.20.Nr,  71.15.Ap,  71.15.Mb 

I.  Introduction  and  Motivations 

Nowotny  et  al  [1]  and  Thomas  [2]  provided  some  of  the  earliest  experimental  data  on 
wurtzite  zinc  oxide  (w-ZnO),  the  former  on  the  lattice  parameters  and  the  latter  on  the 
band  gap.  A  representative  set  of  measurements  of  the  band  gap  of  bulk  w-ZnO  shows  a 
robust  agreement  on  its  low  temperature  value  of  3.44  ±  0.06  eV.  [1-7]  This  range,  i.e., 
±0.06  eV,  is  much  larger  than  experimental  uncertainties;  it  simply  reflects  differences 
that  a  traceable  to  those  of  actual  lattice  constants,  measurement  techniques  and 
temperatures,  and  other  factors.  Some  of  these  factors  (impurities,  defects)  are  associated 
with  the  growth  conditions  of  the  samples.  Srikant  and  Clarke  [8]  appear  to  have  settled 
the  debate  on  the  room  temperature  band  gap  of  w-ZnO.  They  utilized  different 
measurement  techniques  to  arrive  at  a  value  of  3.30  eV.  These  techniques  included 
reflection  and  transmission  absorption,  ellipsometric  spectroscopy,  Fourier  transform 
infrared  (FTIR)  spectroscopy,  and  photoluminescence.  This  multitude  of  measurement 
techniques  avoided  limitations  reported  to  afflict  some  previous  investigations  that 
reported  a  room  temperature  band  gap  of  bulk  w-ZnO  0.1  and  0.2  eV  below  the  currently 
accepted  value  of  3.30  eV.  Other  groups  reported  photoemission  [9-11]  results  for  the 
band  structure,  the  electrons  effective  mass  [12]  and  the  bulk  modulus  [13]. 

Several  groups  [14-20]  have  reported  measured  band  gaps  for  films  of  w-ZnO  grown  by  a 
variety  of  techniques.  These  growth  methods  included  spray  pyrolisis  [14],  pulsed  laser 
deposition  [16,18],  cathodic  electrodeposition  in  aqueous  solutions  [17],  ion  layer  gas 
reaction  (ILGAR)  and  rf  magnetron  sputtering  [19],  and  plasma  enhanced  chemical  vapor 
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LDA)  effective  masses  of  0.25  (0.21),  0.25  (0.21),  and  0.21  (0.18)  mo  in  the  kx,  ky,  and  kz 
directions  respectively,  where  the  values  between  parentheses  resulted  from  their  LDA 
calculations.  These  results  were  obtained  for  experimental  lattice  constants  of  3.24961  A 
,  5.20653  A  and  u  parameter  of  0.345.  Lambrecht  et  al.  [33]  noted  that  their  calculated 
values  of  0.23  mo  and  0.21  mo,  in  the  parallel  and  perpendicular  directions,  are  likely  to 
be  underestimates.  Xu  and  Ching  [35]  reported  0.37,  0.28,  and  0.32  mo  in  the  T-K,  F-A, 
and  F-M  directions  in  the  Brillouin  zone,  for  lattice  constants  given  in  Table  1.  These 
authors  utilized  the  orthogonalized  LCAO  formalism  and  employed  a  local  density 
potential  plus  additional  exchange  derived  with  Wigner’s  interpolation.  The  1966 
measurements  by  Baer  [12]  found  the  electron  effective  mass  for  w-ZnO  to  be  0.24 
±  0.02  mo.  The  above  LDA  results  for  the  electron  effective  mass  are  uniformly  lower 
than  this  experimental  value. 

The  above  overview  of  previous  theoretical  results,  with  the  general  disagreement 
between  calculations,  on  the  one  hand,  and  between  calculations  and  the  firmly 
established  experimental  data,  on  the  other  hand,  provides  the  motivation  for  this  work. 
The  wide  spread  of  the  LDA  results  for  the  band  gap  of  w-ZnO,  from  0.23  to  2.26  eV, 
suggests  that  limitations  other  than  those  genuinely  attributable  to  the  potential  affected 
these  results.  Since  1998,  Bagayoko,  Zhao,  and  Williams  have  introduced  a 
computational  method  [36-39]  that  clearly  showed  that  LDA  potentials  are  no  longer 
known  to  lead  to  woeful  underestimates  of  the  band  gap  of  semiconductors.  Indeed, 
using  the  Rayleigh  theorem,  these  authors  introduced  a  new  form  of  convergence  for 
variational  computations  that  utilize  a  basis  set.  The  need  to  seek  methodically  [36-37]  a 
basis  set  whose  size  is  converged  [i.e.,  the  optimal  basis  set]  vis  a  vis  the  description  of 
the  occupied  states,  they  argued,  straightforwardly  arises  from  (a)  the  use  of  only  the 
wave  functions  of  the  occupied  states  to  construct  the  charge  density,  the  potential,  and 
the  Hamiltonian  in  the  iterative  process  and  (b)  the  Rayleigh  theorem  that  asserts  the 
possible,  continuing  lowering  of  some  unoccupied  energies  or  bands  when  basis  sets 
larger  than  the  optimal  one  are  utilized.  This  extra-lowering  of  some  unoccupied 
energies,  even  though  the  physics  (i.e.,  the  Hamiltonian)  is  no  longer  changing,  was 
identified  as  a  basis  set  and  variational  effect  that  has  afflicted  most  of  the  previous 
calculations,  including  those  utilizing  the  generalized  gradient  approximation  (GGA)  of 
density  functional  theory  (DFT)  [40-41]  and  the  GWA  approximation  [42].  Hence,  the 
clear  motivation  of  this  work  is  to  utilize  properly,  as  per  the  BZW  method,  a  local 
density  potential  to  describe  the  electronic  and  related  properties  of  w-ZnO.  In  the 
remaining  of  this  article,  we  briefly  present  our  method  and  discuss  our  results,  which  are 
basically  in  agreement  with  experiment,  particularly  as  compared  to  previous,  theoretical 
ones. 


n.  Method 

Our  computational  method  has  been  extensively  described  in  previous  publications  [36- 
39].  It  is  characterized  by  the  use  of  a  local  density  approximation  (LDA)  potential  [43] 
as  parameterized  by  Vosko  et  al  [44],  the  well-known  linear  combination  of  Gaussian 
orbitals  (LOGO),  and  the  rigorous  application  of  the  BZW  method  in  carrying  out  the 
calculations.  We  utilized  a  program  package  developed  and  refined  over  decades  [45-46]. 
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III.  Results 


Some  of  our  calculated  band  structures  are  shown  in  Figures  1  and  2.  Figure  1  illustrates 
the  dramatic  differences  between  our  results  and  those  or  previous  LDA  and  other 
calculations.  Essentially,  there  are  no  major  differences  between  our  calculated  occupied 
bands  and  some  previously  reported  ones.  For  the  unoccupied  bands,  however,  our  results 
are  distinctively  far  those  of  previous  calculations.  In  Figure  1,  the  lowest  unoccupied 
energy  at  the  F  point  is  shifted  downward  by  approximately  2.2  eV  as  compared  to  its 
value  obtained  with  the  optimal  basis  set.  Clearly,  the  differences  between  the  low 
energy  conduction  bands  from  calculation  11  and  III  are  far  from  those  expected  from  a 
rigid  shift.  Hence,  they  cannot  be  corrected  with  simplistic  "scissors"  operations. 

The  above  differences  between  our  calculated  conduction  bands  and  those  in  previous 
reports  explain  the  reason  our  calculated  band  gaps  of  3.2  and  3.278  eV,  at  lattice 
constants  specified  in  Table  1  and  Figures  1  and  2,  are  very  close  to  experimentally 
measured  values.  The  differences  between  our  results  from  the  calculation  at  the 
experimental  lattice  constants  in  Figure  2  and  the  one  just  below  it,  in  Table  1,  are  mainly 
due  to  the  difference  in  the  value  of  u,  given  that  the  lattice  parameters  are  not  far  apart. 
Schrder  and  coworkers  found  a  similar  situation  in  their  LMTO  and  GWA  calculations. 
[27]  The  value  of  u,  understandably,  strongly  influences  that  of  the  band  gap. 

Our  calculated,  total  (DOS)  and  partial  density  of  states  (pDOS)  in  Figures  3  and  4 
naturally  reflect  the  differences  discussed  above  for  the  electronic  energy  bands.  In 
particular,  the  location  of  the  lowest  conduction  bands  with  respect  to  the  Fermi  energy  in 
these  figures  is  drastically  different  from  that  in  most  of  the  previous  theoretical 
investigations.  The  inset  in  Figure  3  illustrates  our  concept  of  a  practical  band  gap.  [37] 
Essentially,  due  to  instrumental  sensitivity,  analysis  techniques,  and  related  uncertainties, 
an  experimental  work  could  find  values  of  the  band  gap  from  the  theoretical  minimum  of 
3.278  to  close  to  4.0  eV.  This  possibility,  practically  non-existent  in  highly  sensitive 
photoluminescence  measurements,  increases  to  a  likelihood  if  simple  optical  absorption 
is  the  sole  measurement  technique.  The  determination  of  the  onset  of  the  band  edge 
absorption  is  a  source  of  large  uncertainties,  in  addition  to  the  fitting  often  involved  in  the 
determination  of  the  band  gap. 

The  electron  effective  masses  from  our  calculations  exhibit  a  slight  anisotropy,  as 
expected  in  a  wurtzite  structure.  The  effective  mass  is  a  measure  of  the  curvature  of  the 
calculated  bands.  The  agreement  between  calculated  and  measured  effective  masses 
indicates  an  accurate  determination  of  the  shape  of  the  bands.  We  calculated  the 
effective  masses  of  the  n-type  carriers  of  ZnO,  using  the  electronic  structure  from 
calculation  11  (the  solid  line  in  Fig.  1).  Near  the  bottom  of  the  lowest  conduction  band  at 
the  gamma  point,  we  obtained  0.254,  0.260,  and  0.264  mo  in  the  F-A,  F-K,  and  F-M 
directions,  respectively,  for  the  lattice  constants  in  Figure  1.  For  the  conduction  band  in 
Figure  2,  the  corresponding  values  are  0.257,  0.258,  and  0.257  mo  respectively.  These 
values  agree  very  well  with  the  measured  value  of  0.24  mo,  within  the  experimental 
uncertainty  [12]  of  0.02  mo,  and  with  the  experimental  value  of  0.275  mo.  [48]  The  above 
excellent  agreement  with  experiment  and  the  underestimates  from  previous  LDA 
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Table  1 .  Calculated,  fundamental  band  gaps  (Eg,  in  eV)  of  w-ZnO,  along  with 
pertinent  lattice  constants  in  Angstroms,  compared  to  experiment.  Numbers  in 
the  last  four  columns,  for  a  given  row,  are  from  the  reference  cited  in  that  row. 


Local  Density 

Approximation 

(LDA) 

Potentials 

I  Computational 

Method 

wm 

m 

u 

Eg  (eV) 

LCAO-BZW 
(Present  work) 

3.2501 

5.2071 

0.3817 

3.20 

3.2496 

5.206 

0.345 

3.28 

3.2530 

5.2130 

0.3817 

3.21 

3.2700 

5.1800 

0.381 

3.22 

Zn^"'"  Pseudopotential 
Zn^^'^  Pseudopotential 

3.25 

5.21 

3.25 

5.21 

lEl^^l 

Pseudopotential 

3.23 

5.18 

APW 

LAPW 

3.253 

5.2129 

0.3817 

0.77“ 

FLAPW 

0.93“ 

LMTO 

3.2427 

5.1948 

0.3826 

msmmi 

LMTO 

3.2496 

5.2065 

0.345 

lEEEHH 

LMTO 

3.253 

5.2129 

0.3825 

LMTO  (cubic,  a=4.57  A) 

2.269 

FP-LMTO 

0.382 

HOHI 

LDA  &  Additional 
Exchange 

Orthogonalized  LCAO 

3.249 

5.207 

0.345 

0.88' 

LDA+SIC 

Pseudopotential 

3.29 

5.29 

3.77“ 

X“  Potential 

Pseudopotential 

3.25 

5.21 

■■11 

GWA 

LMTO 

3.2427 

5.1948 

0.3826  1 

4.28  and 
4.06“ 

GWA 

LMTO 

3.2496 

5.2065 

0.345 

3.45  and 
3.63“ 

GWA 

LMTO 

3.253 

5.2129 

0.3825 

GWA 

LAPW 

3.253 

5.2129 

0.3817 

Model  GWA 

FLAPW 

■■11 

Empirical  Pseudopotential  (EMP) 

3.5' 

Empirical  KKR 

3.30'“ 

Experiment  (Bulk,  Low  Temperatures) 

3.44 

Experiment  (Bulk,  Room  Temperature) 

3.30 

^Reference  27  “Reference  28  “Reference  29 

“‘Reference  30  “Reference  23  ^Reference  31 

^Reference  32  ‘^Reference  33  'Reference  35 

^Reference  25  ““Reference  26 
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Figure  2.  Calculated  band  structure  of  wurtzite  zinc  oxide  (w-ZnO),  as  obtained  with 
BZW  method,  for  experimental  lattice  constants  of  a  =  3.2496  A,  c  =  5.206  A,  and 
u  =  0.345  .  The  band  gap  is  3.278  eV. 
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Figure  4.  Calculated  partial  density  of  states  (pDOS)  for  w-ZnO  at  the  experimental 
lattice  constants  of  a  =  3.2496  A,  c  =  5.206  A,  and  u  =  0.345.  The  graph  illustrates  the 
influence  of  the  zinc  d  and  the  oxygen  p  in  the  upper  valence  bands.  Zinc  s  dominates  in 
the  lowest  conduction  band. 
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TABLE  I.  Experimental  and  theoretical  lattice  constants  (a,  in  angstroms)  for  c-InN,  along  with  the  calculated  values  of  the  bulk  modulus  gigapascal  and  of 
the  fundamental  band  gap  (eV),  Results  in  the  last  three  columns,  for  a  given  row',  are  from  the  reference  cited  in  that  row. 


Computational  Method 

a{k) 

B(GPa) 

Es(eV) 

t. 

LCAO-BZW  (Present  work) 

5.017 

4.98 

145 

+0.65 

+0.74 

4.95 

145 

-0.36” 

5.004  , 

140 

-0.40' 

Pseudopotential  Method  (PP) 

4.95 

4.91'‘ 

145' 

4.932 

140 

-0.35' 

Local  Density  Approximation  (LDA)  Potentials 

4.788 

155‘ 

-0.18' 

LAPW 

4.94 

145' 

Full  potential  LAPW 

5.03 

138 

1  1 

o  o 

oo  •— 

Full  Potential  LMTO 

4.92 

1391 

-0.4‘ 

Atomic  sphere  approximation  (ASA) 

-0,l‘ 

+0.02  and  +0.08*’ 

Generalized  gradient  approximation  (GGA) 

PP 

5.06 

5.109 

120' 

118 

-0.55' 

LDA  plus  self-interaction  correction  (SIC) 

5.05' 

+0.43” 

QP  Calculation 

PP 

+0.52“ 

QP+SIC 

PP 

+1.31” 

DFT  Exact  Exchange 

+1.4® 

DPT,  SX 

ASA 

+1.3‘ 

Estimate  of  the  bulk  modulus  of  zinc-blende  indium  nitride  (c-InN)  using  elastic  properties  of  wurtzite  InN 

137‘ 

Empirical  Pseudopotential  Calculations  (EMP) 
Experimental:  Measured  lattice  constants 

4.97±0.01" 

+0.592’” 

4.98' 

4.986'’ 

“Reference  10. 

‘Reference  16. 

'’Reference  13. 

^Reference  17. 

'Reference  11. 

‘‘Reference  19. 

'Reference  15. 

‘Reference  20. 

'Reference  18. 

“Reference  22. 

‘Reference  21. 

"Reference  7. 

^Reference  14. 

"Reference  8. 

'’Reference  12. 

‘Reference  9. 

II.  METHOD  AND  COMPUTATIONAL  DETAILS 

We  performed  zero  temperature,  nonrelativistic  calcula¬ 
tions  of  the  electronic  and  related  properties  of  c-InN.  Our 
ab  inito,  self-consistent  calculations  employed  the  local  den¬ 
sity  approximation  potential  of  Ceperley  and  Alder^"^  as  pa¬ 
rametrized  by  Perdew  and  Zunger.^^  As  stated  above,  we 
used  the  LCAO.  The  feature  distinguishing  our  computa¬ 
tional  method  from  the  previous  investigations  noted  above 
consists  of  our  implementation  of  the  BZW  procedure.  In  so 
doing,  we  started  the  calculations  for  c-InN  with  a  minimal 
basis  set.  We  subsequently  performed  several  other  self- 
consistent  calculations  with  larger  and  larger  basis  sets.  The 
basis  set  for  any  of  these  calculations  was  obtained  by  aug¬ 
menting  the  one  for  the  previous  calculation  with  the  orbital 
describing  the  next  excited  level  of  the  atomic  or  ionic  spe¬ 
cies  present  in  the  system.  The  occupied  bands  of  a  given 
calculation  are  compared  to  those  of  the  previous.  These 
comparisons,  particularly  for  the  first  two  calculations,  often 
show  differences  (in  numerical  values,  branching,  or  curva¬ 
ture).  This  process  continues  until  the  occupied  energies 


from  a  calculation  are  equal,  within  computational  uncertain¬ 
ties,  to  their  corresponding  ones  from  the  calculation  that 
follows  it.  Then,  the  output  of  the  former  calculation  pro¬ 
vides  the  physical  description  of  the  material  under  study  and 
the  related  basis  set  is  dubbed  the  optimal  basis  set.  Accord¬ 
ing  to  the  Rayleigh  theorem, some  of  the  unoccupied  bands 
from  the  latter  may  be  lower  than  their  counterpart  from  the 
former.'^-^®’^®''° 

The  above  additional  lowering  is  the  basis  set  and  varia¬ 
tional  effect  inherently  associated  with  variational  calcula¬ 
tions  of  the  Rayleigh-Ritz  type.  In  the  iterative  process,  the 
use  of  the  wave  functions  for  the  occupied  states  only  in  the 
construction  of  the  charge  density,  and  hence  the  potential 
and  the  Hamiltonian,  ensures  the  exhaustion  of  the  account¬ 
ing  for  the  physicaTinteractions  when  the  occupied  energies 
converge  vis-a-vis  the  size  of  the  basis  set.  As  fully  ex- 
plained  elsewhere,'  “  however,  some  unoccupied  ener¬ 
gies  will  continue  to  be  lowered  as  the  basis  set  is  increased 
beyond  the  optimal  one.  The  sizes  of  the  minimal  and  opti¬ 
mal  basis  sets  vary  vastly  with  the  type  of  functions  in  the 
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Energy  (eV) 

FIG,  3.  pDOS  for  zinc-blende  indium  nitride  (c-InN)  as  obtained  with  the 
bands  shown  in  Fig.  1.  The  dominance  of  nitrogen  p  at  the  top  of  the 
valence  band  is  obvious  in  this  graph. 

the  latest  experimental  one  of  4.986  A.  Our  calculated  bulk 
modulus  of  145  GPa  is  close  to  the  results  from  some  other 
LDA  calculations*^’^^  and  disagrees  with  findings  around 
120  GPa  from  GGA  calculation^^’^^  as  shown  in  Table  I.  The 
tables  for  the  calculated  energies  are  expected  to  provide 
useful  comparisons  for  future  experimental  investigations. 
Our  calculated  LDA-BZW  band  structure  is  drastically  dif¬ 
ferent  from  the-  findings  of  most  of  the  previous 
calculations — as  far  as  the  unoccupied  bands  are  concerned. 
The  large  differences  between  our  calculated  band  gaps  and 
previous  LDA  results  in  Table  I  are  a  direct  consequence  of 
the  differences  between  the  respective  unoccupied  bands. 

The  empirical  pseudopotential  (EMP)  result^^  of 
0.592  eV  is  relatively  cloie  to  our  findings  of  0.65  and 
0.74  eV.  We  cannot  draw  much  fundamental  significance 
from  this  fact,  however,  given  that  this  EMP  result  was  ob¬ 
tained  using  model  potential  parameters  whose  derivation 
entailed  fitting  to  data  that  included  a  band  gap  of  0.59  eV. 
For  the  purposes  of  application,  however,  this  closeness  por¬ 
tends  much  importance.  Indeed,  we  expect  potential  param¬ 
eters  derived  from  fitting  to  our  data  to  lend  themselves  to 
credible  and  practical  descriptions  of  electronic,  optical,  elas¬ 
tic,  and  structural  properties  of  materials.  This  assertion  is 
partly  supported  by  the  versatility  and  relative  ease  of  em¬ 
pirical  pseudopotential  calculations. 

As  in  the  case  of  w-InN,^^  our  results  do  not  show  any 
indication  of  an  overestimation  of  the  p-d  repulsion  by  LDA 
potentials.  This  overestimation  was  believed^”  to  be  the 
source  of  the  very  small  or  negative  band  gaps  by  pushing 
the  top  of  the  valence  band,  dominated  by  p  states,  to  higher 
energies.  According  to  our  findings,  it  is  rather  the  extra  low¬ 
ering  of  the  bottom  of  the  conduction  band  that  produces 
LDA  band  gaps  that  are  negative  or  very  small  ’  if  LCAO 
type  computations  do  not  utilize  the  BZW  approach  to  avoid 
it  while  ensuring  the  adequacy  (or  convergence)  of  the  basis 
set  for  the  description'  of  occupied  states.  In  fact,  the  basic 
derivation  of  the  ground  state  theory  that  is  the  original  den¬ 
sity  functional  theory^'*’^^  implicitly  suggests  such  an  ap¬ 
proach,  notwithstanding  the  need  to  account  for  the  redistri¬ 
bution  of  electrons  in  molecular  to  solid  environments  with 
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judicious  polarization  and  diffuse  orbitals. 


Lattice  Constant  (Angstrom) 

FIG.  4.  Total  energy  of  (c-InN)  vs  the  lattice  constant.  The  total  energy  at 
the  equilibrium  lattice  constant  of  5.017  A  is  1949.4215  eV. 

As  per  their  definition,  effective  masses  provide  a  mea¬ 
sure  of  the  quality  (i.e.,  curvature)  of  band  structures.  Our 
LDA-BZW  calculations  found  electron  effective  masses,  at 
the  bottom  of  the  conduction  band,  of  0.065mo,  0.066mo,  and 
0.066mo  in  the  F-L,  F-X,  and  F-K  directions,  respectively, 
for  the  equilibrium  lattice  constant.  For  a  lattice  constant  of 
4.98  A,  the  corresponding  effective  masses  for  the  electrons 
are  0.07 6mo,  0.073mo,  and  0.073mo  in  the  F-L,  F-X.  and 
F-K  directions,  respectively.  Our  calculated,  equilibrium, 
electron  effective  masses  are  very  close  to  the  0.066mo  from 
the  EMP  calculations  of  Fritsch,  Schmidt,  and  Grundmann.^^ 
This  agreement  supports  our  comment  above  relative  to  the 
potential  use  of  EMP  calculations  when  the  potential  param¬ 
eters  are  derived  in  part  by  fitting  to  LDA-BZW  results. 

IV.  CONCLUSION 

In  lieu  of  a  conclusion,  we  contend  that  our  ab  initio, 
self-consistent  LDA-BZW  calculations  have  predicted  elec¬ 
tronic,  structural,  and  elastic  properties  of  cubic  InN  (c 
-InN)  in  the  zinc-blende  structure.  It  is  hoped  that  experi¬ 
mental  investigations  will  follow  in  the  near  future.  It 
emerges  from  this  work  that  theoretical  efforts  should  be 
directed  to  the  determination  of  actual  limitations  of  LDA 
and  of  other  approaches  as  opposed  to  echoing  the  chorus 
now  known  to  be  untenable,  as  per  the  physical  interactions, 
and  that  ascribes  to  LDA  a  woeful  underestimation  of  band 
gaps  that  is,  as  per  the  Rayleigh  theorem,  straightforwardly  a 
consequence  of  a  basis  set  and  variational  effect. 
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gap  of  1.5  eV.  The  linear  combination  of  Gaussian  orbital 
(LOGO)  calculations  of  Xu  and  Ching^®  employed  a  LDA 
potential  plus  additional  exchange  effects  obtained  with 
Wigner  interpolation  formula.  They  reported  a  band  gap  of 
■  1.02  eV.  The  LDA-fC  work  of  Wei  et  augmented  the 
LDA  potential  with  atom-dependent  corrections  and  obtained 
a  band  gap  of  0.85  eV  with  an  uncertainty  of  0.1  eV. 

Clearly,  the  above  state  of  theoretical  calculations  is  am¬ 
ply  sufficient  to  motivate  the  present  work.  A  recent  experi¬ 
mental  development^®  exacerbated  the  urgency  of  our  theo¬ 
retical  investigations.  Indeed,  the  concordant  findings  of  the 
second  group  of  experiments,  described  above,  have  just 
been  questioned  by  Shubina  et  al?^  in  2004.  These  authors 
showed  the  possibility  of  explaining  the  spectra  of  this  sec¬ 
ond  group  of  experiments  not  in  terms  of  near-band-edge 
absorption  or  emission,  but  rather  as  manifestations  of 
surface/gap  states  or  Mie  resonances  due  to  In  precipitates  in 
InN.  They  noted  the  poor  thermal  stability  of  InN  and  the 
low  In  vapor  pressure  to  support  their  position.  Specifically, 
for  high-quality  InN  films  grown  at  r>500  °C  by  MBE  or 
MOVPE,  they  utilized  x-ray  diffraction  to  observe  tetragonal 
In  precipitates  in  hexagonal  InN.  Their  very  low-temperature 
(0.35  K)  measurements  by  thermally  detected  optical  ab¬ 
sorption  (TDOA)  indicate  the  possibility  of  having  a  band 
gap  around  1 .4  eV.  The  authors  refrained  from  making  a 
definitive  statement  as  to  the  location  of  the  fundamental 
band  gap,  despite  their  observations,  due  to  the  difficulties  in 
separating  the  intraband  absorption  in  InN  from  Mie  reso¬ 
nances. 

METHOD  AND  COMPUTATIONAL  DETAILS 

Our  calculations  employed  the  Ceperley  and  Alder 
local-density-functional  potential  as  parametrized  by  Perdew 
and  Zunger^^  in  an  ab  initio  LOGO  formalism.  We  also  em- 
ployed  the  GGA  potential  of  Perdew  and  co-workers  which 
goes  beyond  LDA.  The  calculations  are  nonrelativistic  and 
for  zero  temperature.  A  key  distinction  between  this  work 
and  previous  calculations  resides  in  our  utilization  of  the  ab 
initio  Bagayoko,  Zhao,  and  Williams^"*’^^  (BZW)  method  in 
carrying  out  the  self-consistent  computations.  Introduced^"*  in 
1998  and  extensively  employed^^  ^  thereafter,  this  method 
rigorously  avoids  a  basis  set  and  variational  effect  inherently 
associated  with  all  variational  calculations  that  utilize  both 
the  linear  combination  of  atomic  orbital  (LCAO)  and  a 
charge  density  obtained  only  with  the  wave  functions  of  the 
occupied  states.  The  effect,  a  consequence  of  the  Rayleigh 
theorem,  consists  of  a  lowering  of  unoccupied  energy  levels 
or  bands,  which  is  not  ascribable  to  a  physical  interaction, 
after  the  total  convergence  of  the  charge  density,  the  poten¬ 
tial,  and  of  the  occupied  energy  levels  with  respect  to  the  size 
of  the  basis  set.^"*’^^ 

BZW  self-consistent  calculations  begin  with  the  mini¬ 
mum  basis  set,  i.e.,  one  that  is  just  large  enough  to  include  all 
the  orbitals  for  the  occupied  levels  in  the  atomic  species  that 
are  present  in  the  system  under  study.  This  initial  basis  set,  in 
our  case,  is  derived  from  calculations  for  the  atomic  or  ionic 
species  present  in  the  finite  or  infinite  system  under  investi¬ 
gation.  Completely  self-consistent  calculations  with  this  ba¬ 


sis  set  are  generally  not  expected  to  describe  the  system  cor¬ 
rectly.  Indeed,  they  generally  do  not  account  for  electron 
redistribution  in  the  system  as  opposed  to  that  in  neutral 
atoms  or  ionic  species.  A  second  self-consistent  calculation  is 
performed  with  a  basis  set  that  includes  the  minimum  basis 
set  as  augmented  with  the  orbital  representing  the  lowest 
excited  levels  in  the  atomic  (or  ionic  species).  The  occupied 
energy  levels  or  bands  of  calculations  I  and  II  are  compared. 

.  One  generally' finds  differences  between  these  energies  from 
the  two  calculations.  For  solids,  these  differences  are  in  nu¬ 
merical  values  and  the  branching  and  curvatures  of  the 
bands.  A  third  calculation  is  performed  by  augmenting  the 
basis  set  of  calculation  11  as  described  above.  If  a  comparison 
of  the  self-consistent,  occupied  energies  of  calculations  II 
and  III  do  not  show  a  complete  agreement,  within  applicable 
errors,  the  process  of  augmenting  the  basis  set  is  continued. 
Eventually,  calculations  N  and  (A^+1)  will  have  the  same 
occupied  energies.  Further,  the  charge  density  and  the  poten¬ 
tial,  from  these  two  calculations,  are,  respectively,  the 
same.^"*’^^  Upon  the  convergence  of  the  variational  calcula¬ 
tions  with  respect  to  the  size  of  the  basis,  the  results  for 
calculation  N  are  selected  as  the  physical  ones.  One  gener¬ 
ally  finds,  indeed,  that  calculation  (N+l)  produced  some  un¬ 
occupied  energies  that  are  lower  than  the  coiresponding  ones 
in  calculation  M  These  differences  that  vary  vastly  with  the 
value  of  k  in  reciprocal  space  are  manifestations  of  the  afore¬ 
mentioned  basis  set  and  variational  effect.^"*"^®  The  basis  set 
for  calculation  N  is  called  the  optimal  basis  set. 

Other  ordinary  features  of  our  calculations  are  thor¬ 
oughly  described  in  the  literature  in  seminal  papers  that  re¬ 
ported'  on  the  development^®”"**  of  the  program  package  and 
its  recent  utilization.^"*’^®'^®  We  provide  below  computational 
details  germane  to  a  replication  of  our  ab  initio  work.  The 
.'experimental,^  low-temperature  wurtzite  lattice  parameters 
utilized  in  our  work  are  fl=3.544A,  c=5.718A,  and  u 
=0.3790.  Preliminary  calculations  for  the  neutral  In  and  N 
atonns  and  the  solid  w-InN  indicated  that  the  self-consistent 
system  is  approximately  In®-^'''N®'®”.  Consequently, 
we  performed  self-consistent  calculations  for  In'"’'  and  N*" 
to  obtain  the  input,  “atomic”  orbitals.  The  minimal  basis 
set  comprised  atomic  orbitals  representing  In*"^ 
{\s2s2pl>s'ip'idAsApAdSs5p)  and  N*"  (Is2j2p).  We  em¬ 
ployed  even-tempered  Gaussian  functions  in  the  construction 
of  the  atomic  orbitals.  The  s  and  p  orbitals  of  In  were  de¬ 
scribed  with  16  even-tempered  Gaussian  functions  with  the 
respective  minimum  and  maximum  exponents  of  0.14  and 
0.2300  X  10®.  The  d  orbitals  for  In  entail  14  even-tempered 
Gaussian,  the  same  as  for  the  ones  s  and  p  at  the  exclusion  of 
the  two  largest  exponents.  Similarly,  the  j  and  p  orbitals  of  N 
were  constructed  using  13  even-tempered  Gaussian  expo-, 
nents  whose  smallest  and  largest  values  were  0.124  215  and 
0.136  000  5  X  10®.  The  computations  were  generally  com¬ 
pletely  self-consistent  after  60  iterations.  The  convergence, 
for  a  given  self-consistent  calculation,  was  reached  when  the 
potential  did  not  change  by  inore  than  10”^  from  one  itera¬ 
tion  to  the  next.  Upon  reaching  self-consistency,  the  error 
made  in  accounting  for  72  electrons  was  0.004  31. 
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TABLE  I.  Eigenvalues  (eV),  at  some  high  symmetry  points,  for  wurtzite  indium  nitride  (w-InN)  as  obtained 
from  LDA-BZW  calculations  for  a=3.544  A,  c=5.718  A,  and  ii=0.3790.  The  Fermi  energy  of  -0.222  12  eV  is 
set  to  zero  in  the  table. 


A 

L 

M 

H 

K 

r 

-13,0972 

-11.5668 

-11.2329 

-11.3021 

-11.2827 

-11.6988 

-3.25532 

-5.25538 

-5.20231 

-4.89999 

-4.32501 

-5.61294 

-3.25532 

-5.25538 

-4.27059 

-4.89999 

-4.32501 

-0.89662 

-0.46069 

-1.62671 

-3.34893 

-3.34349 

-4.2695 

-0.89662 

-0.46069 

-1.62671 

-2.29475 

-3.34349 

-2.39761 

-0.0049 

-0.46069 

-1.53365 

-1.75651 

-1.23894 

-2.39761 

0 

-0.46069 

-1.53365 

-0.8656 

-1.23894 

-2.14019 

0 

3.366891 

4.570732 

4.748424 

6.125875 

6.740313 

0.878663 

3.367163 

4.570732 

6.175673 

6.125875 

6.740313 

3.662409 

the  absorption,  as  per  their  Fig.  2(c)  visibly  commenced 
around  0.9  eV  and  reached  a  plateau  at  2  eV,  in  accordance 
with  our  concept  of  practical  band  gap. 

The  LDA-BZW  effective  masses  for  the  electron  close  to 
the  bottom  of  the  lowest  conduction  band  at  gamma  show  a 
slight  anisotropy.  The  electron  effective  mass  is  around 
0.088mo  in  the  plane  perpendicular  to  fc,  (from  F  to  M  and  F 
to  K)  and  Q.D^lirig,  in  the  direction  of  k,  (from  F  to  A),  in  the 
immediate  vicinity  of  the  gamma  (F)  point.  The  GGA-BZW 
calculations  led  to  values  of  the  electi'on  effective  mass  very 
similar  to  the  ones  for  LDA.  Specifically,  at  the  bottom  of  the 
lowest  conduction  band,  they  are  0.081,  0.077,  and  OMlrrig 
from  F  to  M,  F  to  K,  and  F  to  A  directions,  respectively. 
These  effective  masses  are  also  in  agreement  with 
experiment,  talcing  into  account  the  quality  and  other  fea¬ 
tures  of  the  samples  (i.e.,  temperature,  pressure,  electron 
concentrations).  Wu  et  al.  obtained  a  value  of  0.07  by  using 
a  band  gap  of  0.7  eV  in  their  analysis.  They  also  provided  a 
plot  of  the  electron  effective  mass  as  a  function  of  electron 
concentration. 

DISCUSSIONS 

For  a  discussion  of  our  results,  we  underscore  the  need 
for  specificity  relative  to  the  lattice  parameters,  sample 
quality and  growth  conditions,  temperature, ^\nd  electron 
concentration^^  when  comparing  results  from  different  ex¬ 
perimental  investigations.  In  particular,  the  significant  effect 
of  electron  concentration  on  the  band  gap^^’'*^"^®  could  ex- 
plain  the  0.75-0.8-eV  band  gap  found  by  Wu  et  al.  whose 
samples  had  a  relatively  low  electron  concentration  of  5 
X  10**  cm"^.  Our  calculated  band  gap  of  0.88  eV  agrees  very 
well  with  the  experimental  result  of  0.89  eV  of  Inishima  et 
for  a  sample  with  an  electron  concentration  of  5 
X  lO'®  cm"^.  As  noted  above  in  connection  with  the  “practi¬ 
cal  band-gap,”  measurement  techniques,  the  analysis  of  data, 
and  related  uncertainties  are  critical.  Specifically,  we  reiter¬ 
ate  that  absorption  techniques  could  lead  to  a  value  of  the 
band  gap  of  w-InN  as  large  'as  2  e'V. 

Some  authors  evoked  a  p-d  repulsion  to  explain  the 
small  (generally  negative)  band  gap  from  previous  LDA  cal¬ 
culations.  LDA  is  reported  to  overestimate  this  repulsion  and 
to  push  upward  the  top  of  the  valence  band.  Our  results 
provide  no  indication  of  such  a  limitation  of  LDA.  Rather, 


they  clearly  show  that  the  extra  lowering  of  the  bottom  of  the 
conduction  band,  on  account  of  the  Rayleigh  theorem, 
straightforwardly  leads  to  the  band-gap  problem  in  LDA  and 
other  variational  calculations  that  employ  a  basis  set.  This 
work  partly  responds  to  the  calls  by  Ghani  et  al.‘^  and  by 
Vurgaftman  and  Meyer''^  for  revisiting  the  issue  of  the  band 
gap  of  w-InN  following  their  extensive  compilations  of  ex¬ 
perimental  and  theoretical  results  on  indium  nitride.  Ghani  et 
al.  discussed  the  formation  of  oxynitrides  (InON),  particu¬ 
larly  in  aged  or  annealed  polycrystalline  samples,  as  a  pos¬ 
sible  source  for  the  increase  of  the  band  gap  to  values  around 
2  e'V.  These  polycrystalline  films  are  reported  to  be  suscep¬ 
tible  to  containing  a  high  density  of  oxygen  atoms  at  their 
grain  boundaries. 


CONCLUSION 

In  conclusion,  we  claim  to  have  resolved  the  controversy 
surrounding  the  band  gap  of  wurtzite  InN.  We  did  so  within 
the  local-density  approximation  (LDA)  with  the  BZW 
method.  Unlike  many  schemes  purporting  to  correct  limita¬ 
tions  of  local-density  approximation  and  that  come  with 
other  approximations  whose  effects  are  at  best  ill  understood, 
the  ab  initio  BZW  method  is  rigorously  anchored  on  a  math¬ 
ematical  theorem.  The  above  points  explain  the  reason  this 
communication  is  hoped  to  end  the  echoing  of  the  incorrect 
chorus  claiming  that  it  is  well  known  that  local-density  ap¬ 
proximation  woefully  underestimates  the  band  gap  of  semi¬ 
conductors  and  insulators.  In  the  above  sense,  and  even 
though  this  communication  is  devoted  to  wurtzite  InN,  the 
method  applies  to  other  calculations  that  employ  basis  sets  in 
a  variational  scheme  of  the  Rayleigh-Ritz  type,  as  is  the  case 
in  many  energy  level  (finite  systerns)  or  band  calculations. 
An  utterly  imp'ortant  point  associated  with  the  above  points 
stems  from  the  fact  that  the  proliferation  of  schemes  aimed  at 
“correcting”  limitations  of  LDA  has  started  to  look  like  that 
of  epicycles  in  the  Ptolemaic  model  of  the  solar  system.  A 
seriously  unfortunate  consequence  of  this  situation  consists 
of  diverting  efforts  from  grasping  the  actual  capabilities  and 
limitations  of  LDA  and  those  of  schemes  aimed  at  going 
beyond  the  local-density  approximation  or  beyond  density- 
functional  theory  altogether. 
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Abstract 

In  this  work,  we  reexamined  some  of  the  mathematical  and  phys¬ 
ical  properties  of  the  ab-initio  LCAO  calculations  for  the  electronic 
structure  of  ZnSe,  Ge,  and  GaAs.  The  utilization  of  non-strongly  min¬ 
imal  systems  in  the  self-consistent  ab-initio  calculations  could  lead  to 
a  non-uniformity  in  approaching  the  solution  in  the  Ritz-process.  We 
performed  test  computations  for  the  electronic  structure  of  hydrogen 
atom.  We  have  proposed  that  an  optimum  basis  set  may  be  needed 
so  that  the  calculated  electron  density  is  converged  and  the  significant 
scattering  of  the  Ritz-coefficients  may  be  avoided.  We  have  applied 
the  new  method  to  the  calculations  of  the  electronic  structure  of  ZnSe, 
Ge,  and  GaAs.  Our  calculated  results  of  the  electronic  properties  agree 
well  with  experimental  data. 


1  Introduction 

Ab-initio  density  functional  calculation  is  a  very  powerful  tool  to  study  the 
physical  properties  of  materials,  including  semiconductors,  metals,  surfaces, 
interfaces,  and  others. [1,  2,  3,  4,  5,  6,  7]  However,  previous  density  functional 
calculations  in  local  density  approximation  (LDA)  for  semiconductors  or  in¬ 
sulators  often  underestimated  the  band-gaps  by  30-50  %.[8,  9, 10, 11]  Closely 
related  to  the  band-gap  problem,  the  calculated  effective  masses  of  electrons 
(n-type  carriers)  and  the  optical  properties  of  semiconductors  -from  previ¬ 
ous  LDA  computations-  also  disagree  with  experimental  results.  There  have 
been  some  theoretical  efforts  intended  to  address  these  problems,  including 
the  calculations  that  utilize  nonlocal,  energy- dependent,  non-Hermitian  self¬ 
energy  operators. [12,  13,  14,  15,  16,  17,  18]  Others  implemented  the  exact 
exchange  potentials  or  included  the  effects  of  core  states  via  the  exchange 
diagram. [19,  20,  21]  Aryasetiawan  and  Gunnarsson  reviewed  several  compu¬ 
tational  methods,  including  the  GW  method,  aimed  at  describing  excited 
state  properties. [17].  Johnson  and  Ashcroft  recently  utilized  simplified  ap¬ 
plications  of  the  GW  method  to  make  scissors-type  corrections  to  the  band 
gaps  of  semiconductors.  [18] 
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(2) 


on  the  wave  functions  of  the  occupied  electron  states,  i.  e,, 

=  E4- 1  p 

ki 

where  /^.  is  the  Fermi  distribution  function.  The  solution  of  the  Kohn-Sham 
equation  has  to  satisfy  a  condition  that 

J  p{r)dv  =  N  (3) 

where  N  is  the  total  number  of  electrons  in  the  system.  In  this  work,  we  im¬ 
plemented  the  formalism  of  the  linear  combination  of  atomic  orbitals  method 
to  solve  the  Kohn-Sham  equation.  The  essentials  of  the  mathematical  prop¬ 
erties  of  the  calculations  are  those  of  the  Ritz-process  of  the  variational  cal¬ 
culations.  We  will  follow  the  notations  in  the  discussion  of  the  variational 
theory  according  to  S.  G.  Mikhlin[23]. 

In  the  LCAO  method, [24,  25,  26,  27,  28,  29,  30,  31,  32]  we  solve  the 
Kohn-Sham  equation  self-consistly.  We  expand  the  eigenfunction  of  the 
effective  single  electron  state  of  the  many-atom  system  as  a  linear  combina¬ 
tion  of  the  atomic  wave  functions,  which  is  a  Ritz  process  as  termed  in  the 
variational  theory  according  to  Mikhlin’s  notations. [23] 

=  E  Cam{ki)cl>am{k,  0  (4) 

am 

where  Cam  are  the  expansion  coefBcient  or  the  Ritz-coefficients;  (pam  (^^r) 
are  Bloch  wave  functions  which  are  expressed  in  terms  of  atomic  wave  func¬ 
tions  as 

Garnik,  ^  Warn ~  ~  Rl)  (5) 

where  Uam  is  the  atomic  wave  function  of  the  or-th  state  of  the  m-th  atom 
at  the  position  fm  and  is  the  coordinate  element  in  the  energy  space  (or 
Hilbert  space).  The  solution  for  the  Ritz-coefficients  Cam  of  the  Kohn- 
Sham  equation  is  the  one  that  the  total  energy  functional  Etotaiip)  reaches 
its  minimum. 

2.2  The  Requirement  of  Strongly  Minimal  Systems  in  the 
Ritz  Process 

We  shall  study  the  countable  system  of  the  coordinate  elements  in  the 
Hilbert  space,  where 

{Un}  =Ui,U2,U3,...,Un,...  (6) 
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the  onset  of  the  computational  instabilities  in  the  LCAO  method  may  not  be 
marked  by  an  occurrence  of  the  negative  eigen-value  of  the  overlap  matrix. 

Mikhlin  illustrated  the  problems  of  the  utilization  of  non-strongly  min¬ 
imal  systems  in  Ref.  [23]  by  several  examples.  Generally,  the  non-strongly 
minimal  systems  will  not  reach  a  solution  in  the  limit  of  n  ->  oo  in  the 
Ritz-process.  For  finite  values  of  n,  some  of  the  non-strongly  minimal  sys¬ 
tems  may  approximate  the  Ritz-coefficients  with  a  tendency  to  stabilize 
non-uniformly  for  increasing  n.  Other  non-strongly  minimal  systems  may 
show  that  the  Ritz-coefiicients  scatter  significantly.  The  coefficients  may 
change  abruptly  from  one  approximation  to  the  next.  One  may  not  know 
the  exact  solution  for  the  problem  and  has  the  difficulty  of  identifying  an 
estimate  for  the  error  in  the  approximate  solutions.  There  is  so  far  no  math¬ 
ematical  theorem  that  can  identify  exactly  the  properties  of  the  instabilities 
in  the  utilization  of  non-strongly  minimal  systems  in  the  Ritz-process.  The 
relevant  problems  may  include  those  such  as  the  condition  (s)  that  the  in¬ 
stabilities  may  occur,  the  behavior  of  the  instabilities,  and  the  method(s)  to 
identify  them. 

The  problems  became  more  complex  in  the  applications  of  the  numerical 
computations  of  the  electronic  structure  of  solid  state  materials,  since  most 
of  the  basis  sets  of  various  implementations  of  ab-initio  density  functional 
calculations  are  not  strongly  minimal.  The  true  plane-wave  basis  set  sat¬ 
isfies  the  condition  to  be  strongly  minimal.  However,  the  requirement  of 
the  computational  accuracy  and  the  number  of  the  true  plane- waves  in  the 
solid  state  calculations  are  usually  beyond  the  limitation  of  our  computer 
resource  at  present.  The  orthogonalized  plane  waves,  which  are  derived  by 
the  Schmitdt  process  to  be  orthogonal  to  the  core  wave  functions,  are  not 
strongly  minimal. 

One  may  propose  to  orthogonalize  the  atomic  orbitals  for  the  solid  state 
calculations.  However,  if  the  basis  set  is  not  strongly  minimal,  the  orthog- 
onalization  of  the  basis  set  in  solid  state  calculations  may  encounter  the 
similar  instability  problem.  The  orthogonalization  procedure  of  the  basis 
set  involves  the  computation  of  the  inverse  of  the  overlap  matrix.  If  the 
smallest  eigen-value  of  the  overlap  matrix  approaches  zero  or  negative,  an 
instability  problem  may  occur. 

2.3  Consideration  of  Some  Physical  Properties 

Although  the  true  plane-wave  basis  set  can  be  implemented  in  ab-initio 
density  functional  calculations  of  the  electronic  structure  of  solid  state  ma- 
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atomic  wave  functions,  we  performed  ab-initio  atomic  calculations  that  also 
utilized  density  functional  formalism.  The  radial  parts  of  the  atomic  wave 
functions  were  expanded  in  terms  of  a  linear  combination  of  Gaussian  or¬ 
bitals  (LOGO). 

The  electronic  energy  levels  and  related  wave  functions  of  hydrogen  atom 
are  exactly  known  from  the  analytical  solution  of  the  Schrddinger  wave  equa¬ 
tion.  We  can  test  the  numerical  performance  of  the  variational  computations 
of  the  Ritz-process  on  hydrogen  atom.  We  utilized  the  same  ab-initio  atomic 
computation  program  and  the  LOGO  formalism  to  calculate  the  electronic 
structure  of  hydrogen  atom.  The  Schrddinger  wave  equation  for  the  calcula¬ 
tion  of  the  electronic  structure  of  hydrogen  atom  is  a  one-electron  problem, 
instead  of  the  many-body  problems  such  as  those  for  the  studies  of  solid 
state  materials.  In  the  ab-initio  computation  of  the  hydrogen  atom  using 
the  LOGO  method,  we  did  not  include  the  density  functional  potentials 
in  the  Hamiltonian.  The  Schrddinger  wave  equation  for  hydrogen  atom  is 
then  a  linear  differential  equation,  instead  of  the  non-linear  differential  equa¬ 
tion  such  as  the  density  functional  Kohn-Sham  equation.  We  compared  the 
results  of  the  numerical  computations  of  hydrogen  atom  with  the  exact  solu¬ 
tion.  In  the  numerical  computations,  we  utilized  basis  sets  of  even-tempered 
Gaussian  exponentials  with  a  minimum  of  0.10  x  10~^  and  a  maximum  of 
0.15  X  10^  in  atomic  unit.  We  used  various  sizes  of  the  basis  sets  of  the 
Gaussian  functions. 

In  Fig.  1,  we  present  the  calculated  electronic  energy  levels  of  hydrogen 
atom  versus  the  total  numbers  of  the  utilized  Gaussian  functions  (Gaussian 
orbitals).  The  solid,  dotted,  and  dash-dotted  lines  in  Fig.  1  present  the 
calculated  energy  levels  of  s-,  p-,  and  d-states  of  hydrogen  atom,  respectively. 
As  the  numbers  of  the  Gaussian  functions  in  the  basis  sets  increase  from  5 
to  15,  the  calculated  energy  levels  approach  their  converged  values.  The 
calculated  first  three  energy  levels  of  hydrogen  atom  converged  to  -1.0,  - 
0.250,  and  -0.111  Ry,  when  we  utilized  15  -  60  Gaussian  functions  in  the 
basis  sets.  These  calculated  energy  levels  of  hydrogen  atom  agree  very  well 
with  that  of  the  exact  solution  for  the  principal  quantum  number  n  =  1,  2, 
and  3,  respectively.  The  lowest  energy  level  in  Fig.  1  belongs  to  the  Is  state 
of  H  atom.  The  second  lowest  level  includes  2s  and  2p  states.  The  third  one 
includes  3s,  3p,  and  3d  states. 

When  the  numbers  of  the  utilized  Gaussian  functions  in  the  basis  sets 
are  larger  than  about  70,  the  calculated  electron  energy  levels  of  the  3p  and 
3d  states  of  hydrogen  atom  dropped  to  lower  values  than  those  of  the  exact 
solution.  The  features  illustrate  the  instabilities  of  the  variational  calcu- 
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of  the  five  self-consistent  calculations.  We  first  carried  out  completely  self- 
consistent  calculations  for  ZnSe  using  a  basis  set  that  included  Zn(ls2s3s4s 
2p3p  3d)  and  Se(ls2s3s4s  2p3p4p  3d)  orbitals.  We  then  repeated  the  self- 
consistent  calculation  using  the  above  basis  set  as  augmented  by  the  orbitals 
describing  the  next  lowest-lying  energy  level  of  Zinc.  Hence,  Zn(4p)  orbitals 
were  added  to  basis  set  I.  We  then  plotted  the  energy  bands  from  these 
two  distinct  calculations,  i.e.,  calculation  I  and  II,  and  compared  them  nu¬ 
merically  and  graphically.  Differences  were  obvious.  We  then  performed 
calculation  III  and  compared  its  results  to  those  of  calculation  II.  Up  to  five 
self-consistent  calculations  were  performed  and  the  results  of  each  calcula¬ 
tion  were  compared  to  those  of  the  previous  one  as  explained  above. 

In  Fig.  2,  we  presnt  the  calculated  electron  energy  bands  of  ZnSe  from 
calculation  III  (solid  lines)  and  calculation  IV  (dashed  lines).  Here  F  = 
(0, 0, 0);  D  =  (1, 1,  l)7r/a;  ^  =  (0, 1, 0)27r/a;  K  =  (3/4,  3/4, 0)27r/a.[35]  The 
zero  of  the  energy  was  set  at  the  top  of  the  valence  bands.  The  calculated 
energies  of  the  occupied  valence  bands  and  the  unoccupied  conduction  bands 
converged  to  an  average  difference  of  several  meV.  The  fourth  calculation 
(the  dashed  line)  of  Fig.  2  gives  sufficiently  converged  electron  energy  bands 
of  ZnSe  with  respect  to  the  size  of  the  basis  set.  As  we  present  in  the  follow¬ 
ing  sections,  the  calculated  electronic  properties  of  ZnSe  from  calculation 
IV  agree  well  with  experimental  data. 

However,  as  we  added  more  atomic  orbitals  in  the  calculations,  the  cal¬ 
culated  electron  energy  bands  of  ZnSe  drifted  away  from  the  converged  re¬ 
sults  of  calculation  IV  around  the  F-point.  One  may  observe  a  computa¬ 
tional  behavior  that  is  similar  to  that  of  the  calculations  for  hydrogen  atom. 
Fig.  3  shows  the  comparison  of  results  from  calculation  IV  (solid  lines)  and 
from  calculation  V  (dashed  lines).  The  calculated  occupied  energy  levels 
converged  with  respect  to  the  size  of  the  basis  set.  However,  some  of  the 
calculated  unoccupied  energy  levels  become  lower  than  those  of  the  calcu¬ 
lation  IV,  as  more  atomic  orbitals  are  added  in  the  calculations.  Different 
from  the  case  of  hydrogen  atom,  there  was  no  an  exact  solution  for  the 
electronic  structure  of  ZnSe  as  a  reference.  Various  physical  interpretations 
were  devised  to  explain  the  numerical  performance  of  the  Ritz-process  in 
the  calculations.  Some  of  them  were  believed  to  contribute  in  various  as¬ 
pects  to  the  problem.  In  fact,  the  problem  was  previously  claimed  as  the 
failure  of  density  functional  theory  or  the  failure  of  local  density  functional 
calculations  in  many  occasions.  Part  of  the  support  for  this  claim  was  due 
to  the  fact  that  the  calculated  conduction  band  energies  (or  the  band-gap 
energies)  of  semiconductors  utilizing  various  computational  methods  often 
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of  states  (DOS)  of  ZnSe  in  Fig.  4,  as  obtained  from  calculation  IV.  The 
total  DOS  curve  in  Fig.  4,  particularly  the  inset,  indicates  that  the  small 
values  of  the  density  of  states  from  2.6  eV  to  2.9  eV  strongly  suggest  that 
measurements  or  related  analysis  that  are  not  very  sensitive  may  not  detect 
the  smallest  energy  gap. 

Our  calculated  effective  mass  of  n-type  carriers,  m* ,  near  the  F-point  is 
0.15  mo,  where  mo  is  the  free-electron  mass.  This  result  is  in  good  agreement 
with  the  experimental  value  of  0.16mo.[40,  39]  The  calculated  effective  mass 
of  n-type  carriers,  away  from  the  F-point,  is  between  0.15mo  and  0.17mo. 

We  also  performed  the  self-consistent  ab-initio  calculations  for  the  elec¬ 
tronic  structure  of  ZnSe,  utilizing  local  density  approximation  for  the  exchange- 
correlation  potential.  We  compared  the  electron  energy  bands  from  the  LDA 
and  GGA  calculations.  The  average  difference  between  the  resulted  bands 
of  the  two  calculations  is  of  the  order  of  1  mRy. 

3.3  Electronic  Structure  of  Ge 

We  performed  the  self-consistent  ab-initio  calculations  for  the  electronic 
structure  of  crystalline  Ge  in  the  diamond  structure.  We  used  a  zero  temper¬ 
ature  lattice  constant  of  a  =  5.63A  in  these  calculations.  [39]  We  tested  the 
convergency  of  the  electronic  structure  calculation  with  respect  to  the  size 
of  the  basis  sets  of  atomic  orbitals.  We  calculated  the  atomic  wave  function 
using  a  separate  computer  program  that  employed  ab-initio  density  func¬ 
tional  computations  for  atoms.  We  expanded  the  atomic  wave  functions  as 
a  linear  combination  of  Gaussian  orbitals  (LOGO)  in  real  space  as  discussed 
in  the  previous  section. 

In  the  calculation  of  the  electronic  structure  of  Ge,  we  employed  non¬ 
local  density  functional  potentials  in  the  generalized  gradient  approximation 
that  was  developed  by  Perdew  and  Wang. [22]  We  first  carried  out  com¬ 
pletely  self-consistent  calculations  for  crystalline  Ge  using  a  basis  set  that 
included  atomic  orbitals  of  Ge(ls2s3s4s  2p3p4p  3d).  We  then  repeated  the 
self-consistent  calculation  using  an  augmented  basis  set  that  also  included 
the  orbital  of  Ge(4d°).  In  the  self-consistent  calculations  of  the  electronic 
structure  of  solid  Ge,  we  considered  Ge(ls2s  2p)  as  the  core  states.  All  other 
states  were  treated  as  valence  states  and  were  allowed  to  relax  in  the  self- 
consistent  calculations.  We  then  plotted  the  resulting  energy  bands  from 
these  two  distinct  self-consistent  calculations.  We  observed  that  the  occu¬ 
pied  and  unoccupied  bands  from  the  two  calculations  differ  considerably. 

The  next  step  was  to  repeat  our  procedure  for  the  third,  fourth,  and 
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3.4  Electronic  Structure  of  GaAs 

We  also  applied  the  same  calculation  procedure  as  discussed  in  the  previ¬ 
ous  sections  to  the  study  of  the  electronic  structure  of  GaAs  in  the  zinc 
blende  structure.  We  used  the  experimental  lattice  constant  of  a  =  5.65A 
in  the  calculations. [43,  44]  We  calculated  the  atomic  wave  functions  of  Ga 
and  As  self-consistently.  We  evaluated  the  charge  transfer  in  GaAs  from  the 
self-consistent  ab-initio  calculations  of  the  electronic  structure.  Ga  atoms 
in  GaAs  lose  about  0.3  to  0.4  electrons  per  atom  to  As  atoms.  In  the  self- 
consistent  ab-initio  calculations,  we  employed  non-local  density  functional 
potentials  in  the  generalized  gradient  approximation.  We  used  the  Is,  2s, 
and  2p  states  of  Ga  and  As  as  the  core  states  in  the  frozen  core  approxi¬ 
mation,  and  allowed  all  other  states  to  relax  in  the  self-consistent  ab-initio 
calculations.  In  Fig.  6,  we  present  the  calculated  electron  energy  bands 
of  GaAs  along  some  high  symmetry  directions  in  the  Brillouin  zone.  Here 
r  =  (0,0,0);  L  =  (l,l,l)7r/a;  X  =  (0, 1, 0)27r/a;  K  =  (3/4, 3/4,  0)27r/a. 
The  zero  of  the  energy  was  set  at  the  top  of  the  valence  bands.  In  Fig.  6, 
the  solid  lines  represent  the  calculated  results  using  the  atomic  orbitals  of 
Ga(ls2s3s4s  2p3p4p  3d4d)  and  As(ls2s3s4s  2p3p4p  3d4d),  where  Ga(4d) 
and  As(4d)  are  the  unoccupied  atomic  orbitals.  The  dashed  lines  represent 
the  calculated  bands  for  the  basis  set  that  has  an  additional  atomic  orbital 
of  As(5s).  The  calculated  electron  energy  bands  of  the  occupied  states  of 
GaAs  from  the  two  different  calculations  agree  very  well.  The  average  dif¬ 
ference  of  the  energies  of  the  occupied  states  from  these  two  calculations  is 
about  several  meV. 

GaAs  has  a  direct  band- gap  at  the  F-point  as  shown  in  Fig.  6.  The  lowest 
energy  state  of  the  conduction  band  is  at  the  F-point.  The  lowest  energy  of 
the  conduction  band  at  the  L-point  is  higher  than  that  at  the  F-point  by  0.2 
eV.  This  result  agrees  well  with  the  experimental  value  of  0.2  to  0.3  eV.[44] 
The  theoretical  band-gap  energy  of  GaAs  is  1.24  eV  from  this  calculation. 
The  reported  experimental  values  of  the  band-gap  energy  of  GaAs  is  1.4  to 
1.5  eV.[43,  44,  45]  As  we  discussed  in  the  previous  sections,  we  present  below 
the  calculated  electron  density  of  states  of  GaAs  in  Fig.  7.  There  is  a  tail 
structure  near  the  conduction  band  edge.  The  small  values  of  the  density  of 
states  from  1.24  eV  to  1.7  eV  strongly  suggest  that  measurements  or  related 
analyses  that  are  not  very  sensitive  may  not  detect  the  smallest  energy  gap. 

We  also  calculated  the  effective  mass  m*  of  the  electrons  (n-type  carriers) 
in  the  lowest  conduction  band  near  the  F-point.  The  effective  mass  m* 
is  calculated  from  the  curvature  of  the  band  dispersions.  The  calculated 
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Table  2;  Calculated  electronic  energies  (eV)  of  Ge  at  the  selected  high  sym¬ 
metry  points.  The  experimental  data  are  from  Reference  [39,  44]. 
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Calculation 

Deg 
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Tiu 

-12.59 
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-12.6,  -12.9±0.2 
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0.0 
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0.0 

rL 

0.83 
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0.89 
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2.81 
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3.01 

^2v 

-10.41 
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-10.6±0.5 

Liv 

-7.44 

-7.7±0.2 

Ll 

-1.39 

-1.4±0.3 

Lie 

0.62 

1 

0.74 
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4.02 
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4.3±0.2,  4.2T0.1 
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8.29 

1 

7.8T0.6,  7.9±0.1 
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-8.64 
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-9.3±0.2 

Xiv 

-3.00 
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-3.15±0.2,  -3,5±0.2 

Xlc 

1.04 

2 

1.3±0.2 

Xzc 

9.98 

2 

E,  1 

0.62 

1 

0.66,  0.74 

Table  3:  The  effective  masses  (in  mo)  of  the  n-type  carriers  at  the  lowest 
conduction  band  of  Ge  in  the  diamond  structure. 
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Figure  2:  The  calculated  electronic  energy  bands  of  ZnSe  along  the  high 
symmetry  directions.  The  solid  and  dashed  lines  represent  the  results  of 
calculations  III  and  IV,  respectively. 
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Figure  4:  The  calculated  density  of  states  of  ZnSe.  The  inset  shows  the  tail 
structure  of  the  density  of  states  of  ZnSe  near  the  bottom  of  the  conduction 
bands. 
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Figure  6:  The  calculated  electron  energy  bands  of  GaAs. 
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